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Abstract. Recent developments in the physics of extra dimensions have opened up 
new avenues to test such theories. We review cosmological aspects of brane world 
scenarios such as the Randall-Sundrum brane model and two-brane systems with a 
bulk scalar field. We start with the simplest brane world scenario leading to a consistent 
cosmology: a brane embedded in an Anti-de Sitter space-time. We generalise this 
setting to the case with a bulk scalar field and then to two-brane systems. 

We discuss different ways of obtaining a low-energy effective theory for two-brane 
systems, such as the moduli space approximation and the low-energy expansion. A 
comparison between the different methods is given. Cosmological perturbations are 
briefly discussed as well as early universe scenarios such as the cyclic model and the 
born-again brane world model. Finally we also present some physical consequences 
of brane world scenarios on the cosmic microwave background and the variation of 
constants. 
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1. Introduction 

Whilst theories formulated in extra dimensions have been around since the early 
twentieth century, recent developments have opened up new avenues which have enabled 
the cosmological consequences of such theories to be extracted. This, together with 
advances in high precision cosmological data, opens up the exciting possibility of testing 
and constraining extra-dimensional theories for the first time. In this review we explore 
these developments. 

In the early twentieth century Kaluza and Klein [1], following Nordstrom, 
attempted to unify electromagnetism and gravity by proposing a theory in five 
space-time dimensions, with the extra spatial dimension compactified on a circle. 
More recently, attempts to construct a consistent theory incorporating both quantum 
mechanics and gravity have been constructed in more than four dimensions, with the 
extra, spatial dimensions compactified. In particular superstring theories, which include 
both quantum theory and gravity, are only consistent in ten dimensions j2], with the 
extra dimensions being spatial. The usual four dimensional physics is retrieved by 
compactifying on a compact manifold with dimensions typically of the order of the 
Planck scale. This is essentially a generalisation of Kaluza-Klein theory. 

Another approach has emerged recently, motivated by M-theory. M-theory is an 
umbrella theory that, in certain limits, reduces to the five known string theories or to 
supergravity. Horava and Witten [3] showed that the strong coupling limit of the E 8 x E 8 
heterotic string theory at low energy is described by eleven dimensional supergravity 
with the eleventh dimension compactified on an orbifold with Z 2 symmetry, i.e. an 
interval. The two boundaries of space-time (i.e. the orbifold fixed points) are 10- 
dimensional planes, on which gauge theories (with E 8 gauge groups) are confined. Later 
Witten argued that 6 of the 11 dimensions can be consistently compactified on a Calabi- 
Yau threefold and that the size of the Calabi-Yau manifold can be substantially smaller 
than the space between the two boundary branes [3], i.e. Planckian. The branes are 
hypersurfaces embedded in the extra dimension. Thus, in this limit space-time looks 
five-dimensional with four dimensional boundary branes [S] . This has led to the concept 
of brane world models. 

More generally speaking, in brane world models the standard model particles are 
confined on a hypersurface (called a brane) embedded in a higher dimensional space 
(called the bulk). Only gravity and other exotic matter such as the superstring dilaton 
field can propagate in the bulk. Our universe may be such a brane-like object. This 
idea was originally motivated phenomenologically (see jEJ-jHi]) and later revived by 
the recent developments in string theory. Within the brane world scenario, constraints 
on the size of extra dimensions become weaker because the standard model particles 
propagate only in three spatial dimensions. Newton's law of gravity, however, is sensitive 
to the presence of extra-dimensions. Gravity is being tested only on scales larger than 
a tenth of a millimeter and possible deviations below that scale can be envisaged. 

Arkani-Hamed, Dimopoulos and Dvali (ADD) ^T] and ^2|, following an earlier 
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idea by Antoniadis [13], proposed an interesting, but simple model, considering a flat 
geometry in (4 + <i)-dimensions, in which d dimensions are compact with radius R 
(toroidal topology). All standard model particles are confined to the brane, but gravity 
can explore the extra d dimensions. This gives rise to a modification of the gravitational 
force law as seen by an observer on the brane. Two test masses, m\ and ul%, at distances 
r apart will feel a gravitational potential of 

V(r) « \J r <i? 1 

Ml + l r d+1 

fund 

and 

where Mf un d is the fundamental Planck mass in the higher dimensional space-time. This 
allows one to identify the four-dimensional Planck mass Mpi, 

Ml, = M?+jR d . (3) 

^From above we see that on scales larger than R gravity behaves effectively as the 
usual four-dimensional gravity. However, on scales less than R there are deviations 
and gravity looks truly 4 + d dimensional. Since gravity is tested only down to sizes of 
around a millimeter, R could be as large as a fraction of a millimeter. 

The most important result of the ADD proposal is a possible resolution to the 
hierarchy problem, that is the large discrepancy between the Planck scale at 10 19 
GeV and the electroweak scale at 100 GeV. The fundamental Planck mass could be 
comparable to the electroweak scale as long as the volume of the extra dimensional 
space is large enough. To realise this, their proposal requires more than one extra 
dimension. 

However, progress was made by Randall and Sundrum, who considered non-flat, i.e. 
warped bulk geometries [Ti] . [To] . Their models were formulated in one extra dimension 
with the bulk space-time being a slice of Anti-de Sitter space-time, i.e. a space-time 
with a negative cosmological constant. They proposed a two-brane model in which the 
hierarchy problem can be addressed. The large hierarchy is due to the highly curved 
AdS background which implies a large gravitational red-shift between the energy scales 
on the two branes. In this scenario, the standard model particles are confined on a brane 
with negative tension sitting at y = r c , whereas a positive tension brane is located at 
y = 0, where y is the extra spatial dimension. The large hierarchy is generated by the 
appropriate inter-brane distance, i.e. the radion. It can be shown that the Planck mass 
Mpi measured on the negative tension brane is given by (k = y — A5K5/6), 

Mp\ « e 2krc Ml/k, (4) 

where M 5 is the five-dimensional Planck mass and A 5 the (negative) cosmological 
constant in the bulk. Thus, we see that, if M5 is not very far from the electroweak 
scale M\y ^TeV, we need kr c fa 50, in order to generate a large Planck mass on our 
brane. Hence, by tuning the radius r c of the extra dimension to a reasonable value, one 
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can obtain a very large hierarchy between the weak and the Planck scale. Of course, 
a complete realisation of this mechanism requires an explanation for such a value of 
the radion. In other words, the radion needs to be stabilised at a certain value. The 
stabilisation mechanism is not thoroughly understood, though models with a bulk scalar 
field have been proposed and have the required properties jTHj. 

This two brane model was unrealistic because standard model particles are confined 
to the negative tension brane. However, Randall and Sundrum realised that due to the 
curvature of the bulk space time, Newton's law of gravity could be obtained on a positive 
tension brane embedded in an infinite extra dimension with warped geometry. Small 
corrections to Newton's law are generated and constrain the possible scales in the model 
to be smaller than a millimetre. This model does not solve the hierarchy problem, but 
has interesting cosmological implications. 

A spectacular consequence of brane cosmology, which we discuss in detail in a later 
section, is the possible modification of the Friedmann equation at very high energy |17j . 
This effect was first recognised [IB] in the context of inflationary solutions. As we will 
see, Friedmann's equation has, for the Randall-Sundrum model, the form (^U] and [2U] ) 

H 2 = -V H -p + A, (5) 

relating the expansion rate of the brane H to the (brane) matter density p and the 
(effective) cosmological constant A. The cosmological constant can be tuned to zero by 
an appropriate choice of the brane tension and bulk cosmological constant, as in the 
Randall-Sundrum case. Notice that at high energies, for which 
„ 9QttG n 

P> — , (6) 

where k\ is the five dimensional gravitational constant, the Hubble rate becomes 

H oc p, (7) 
while in ordinary cosmology H oc yfp. The latter case is retrieved at low energy, i.e. 

P«— A (8) 

Of course modifications to the Hubble rate can only be significant before nucleosynthesis, 
though they may have drastic consequences on early universe phenomena such as 
inflation. 

A natural extension to the Randall-Sundrum model is to consider a bulk scalar field. 
This occurs automatically in supersymmetric extensions and makes more contact to 
string theory, which has the dilaton field. The addition of a bulk scalar field allows other 
fundamental problems to be addressed, for instance that of the cosmological constant. 
Attemps have been made to invoke an extra-dimensional origin for the apparent (almost) 
vanishing of the cosmological constant. The self-tuning idea [21] advocates that the 
energy density on our brane does not lead to a large curvature of our universe. On the 
contrary, the extra dimension becomes highly curved, preserving a flat Minkowski brane 
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with apparently vanishing cosmological constant. Unfortunately, the simplest realisation 
of this mechanism with a bulk scalar field fails due to the presence of a naked singularity 
in the bulk. This singularity can be shielded by a second brane whose tension has to be 
fine-tuned with the original brane tension [22] • In a sense, the fine tuning problem of 
the cosmological constant reappears through the extra dimensional back-door. However, 
in supersymmetric extensions of the Randall-Sundrum model the bulk scalar field can 
give rise to a late-time acceleration of the universe once supersymmetry is broken on 
the physical brane world. The naked singularity mentioned above is automatically 
beyond the second brane. There is still a fine-tuning problem here, namely that of 
supersymmetry breaking. We discuss this in detail in later sections. 

One might hope that brane world models could explain the origin of the universe 
itself. Attempts have been made in this direction by considering brane collisions. In 
this scenario the universe undergoes a series of big bangs and big crunches as the branes 
go through cycles of moving apart and colliding. This is an exciting development which 
we briefly discuss in the last section. 

In this review we study the physics of one brane systems both with an empty bulk 
and with a bulk scalar field. Birkhoff 's theorem is discussed and examined in the case 
of a bulk scalar field. We then turn our attention to the case of two brane systems. 
The low-energy behaviour of two brane systems is emphasised. Throughout this review 
we discuss the cosmological implications and predictions of brane worlds. In particular, 
the CMB predictions of simple brane models are reviewed. Finally, we conclude by 
considering brane collisions both in the cyclic and born-again brane world scenarios. 
The review is an extension of an earlier review j2H|, with recent developments and new 
material included. There are other reviews on brane worlds and their cosmological 
implications (see e.g. [21] - [H2]), each with a different emphasis. 

2. A brane in an anti— de Sitter bulk space— time 

We begin with a discussion of the simplest non-trivial (i.e. different from Minkowski 
space-time) bulk, the Anti-de Sitter (AdS) space-time. The AdS space-time is 
a maximally symmetric solution to Einstein's equation with negative cosmological 
constant. 

The first well motivated brane world model with an Anti-de Sitter bulk space- 
time was suggested by Randall and Sundrum. Initially, they proposed a two-brane 
scenario in five dimensions with a highly curved bulk geometry as an explanation for 
the large hierarchy between the Planck scale and the electroweak energy scale j!4j . In 
their scenario the standard model particles live on a brane with (constant) negative 
tension, whereas the bulk is a slice of Anti-de Sitter (AdS) space-time. In the bulk 
there is another brane with positive tension. This is the so-called Randall-Sundrum 
I (RSI) model, in order to distinguish this setup from a one brane model, which they 
proposed immediately afterwards. Analysing the solution of Einstein's equation on the 
positive tension brane and sending the negative tension brane to infinity, an observer 
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confined to the positive tension brane recovers Newton's law if the curvature scale of 
the AdS space is smaller than a millimeter JJ]. Note that the higher-dimensional space 
is non-compact, which must be contrasted with the Kaluza-Klein mechanism, where all 
extra-dimensional degrees of freedom are compact. This one-brane model, the so-called 
Randall-Sundrum II (RSII) model, will be our main focus in this section. 

Before discussing the cosmological consequences of the second Randall-Sundrum 
model, let us first look at the static solution, introducing some valuable techniques. The 
(total) action consists of the Einstein-Hilbert action and the brane action, which have 
the form 

R 



Seh = - [d 5 x^g^(^ + A 5 ), (9) 



S br ane= / ^X^^ \-&) . (10) 



Here, A5 is the bulk cosmological constant and a is the (constant) brane tension. K5 is 
the five-dimensional gravitational coupling constant. The brane is located at y — and 
we assume a Z2 symmetry, i.e. we identify y with —y. The ansatz for the metric is 

ds 2 = e- 2K(y) ri^dx»dx u + dy 2 . (11) 
Einstein's equations, derived from the action above, give two independent equations: 

6K' 2 = - n 2 5 A 5 

3K" = 4*5(y). 
The first equation can be easily solved: 




K = K(y) = \-^A 5 y = ky, (12) 



which tells us that A5 must be negative. If we integrate the second equation from — e 
to +e, take the limit e — > and make use of the Z 2 symmetry, we get 

6K'\ = K 2 5 a (13) 

Together with eq. (|T2*j) this tells us that 

A 5 = -^-o 2 (14) 
o 

Thus, there must be a fine tuning between the brane tension and the bulk cosmological 
constant for static solutions to exist. 

Randall and Sundrum have shown that there is a continuum of Kaluza-Klein modes 
for the gravitational field, which is different from the case of a periodic extra dimension 
where a discrete spectrum is predicted. The continuum of Kaluza-Klein modes lead to 
a correction to the force between two static masses on the brane. To be more specific, 
it was shown that the potential energy between two point masses confined on the brane 
is given by 



G^m\m2 / , , I 2 
T2 



v(r)= r - i + ^ + (is) 
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As gravitational experiments show no deviation from Newton's law of gravity on length 
scales larger than a millimeter [33], I has to be smaller than that length scale. 

A remark is in order: we assumed a Z 2 symmetry, i.e. we are considering the 
orbifold Si/Z 2 . Most brane world scenarios assume this symmetry. The physical 
motivation for considering this orbifold compactification comes from heterotic M-theory 
(as mentioned in the introduction). Of course, assuming the Z 2 symmetry leads to a 
significant simplification of the equations and therefore simplicity is a good motivation, 
too. 

In the following we derive Einstein's equation for an observer on the brane. 



2.1. Einstein's equations on the brane 

There are two ways of deriving the cosmological equations and we will describe both 
of them below. The first one is rather simple and makes use of the bulk equations 
only. The second method uses the geometrical relationship between four-dimensional 
and five-dimensional quantities. We begin with the simpler method. 



2.1.1. Friedmann's equation from five-dimensional Einstein equations In the following 
subsection we will set k 5 = 1 as this simplifies the equations a lot. We write the bulk 
metric as follows: 



ds' 



a 2 b 2 (dt 2 - dy 2 ) - a 2 5 ij dx i dx j . (16) 

This metric is consistent with homogeneity and isotropy on the brane located at y = 0. 
The functions a and b are functions of t and y only. Furthermore, we have assumed flat 
spatial sections, but it is straightforward to include a spatial curvature. For this metric, 
Einstein equations in the bulk read: 

) ci 2 ab a" a'b' 
a 2 ab 

„/2 n, 



a 2 b 2 G\ 



+ — + W 
a ab 



2 1.2 



a z b 



PB + pS(y - y b ) 



a 2 b 2 G% 
a 2 b 2 G\ 



ab 
ab 



a!V_ 
ab 



+ kb 2 



-a 2 b 2 Tl 



ab' a'b 
— + 2— + — + — 



aa 



a 



a* 



ab ab 



-a 2 b 2 Tl 



a 2 b 2 G\ = 



a b b 2 a" b" V 
a b tr abb 1 



-a'b 



1i2 



PB+pS{y - yb) 



J- 



(17) 
(18) 
(19) 

(20) 



where the bulk energy-momentum tensor T fc a has been kept general here. For the 



Randall-Sundrum model we will now take ps 



-Pb 



Ak and TP = 0. Later we 



will make use of these equations to derive Friedmann's equation with a bulk scalar field. 
In the equations above, a dot represents the derivative with respect to t and a prime a 
derivative with respect to y. 
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If one integrates the 00-component over y from — e to e and use make use of the Z 2 
symmetry (which implies here that a(y) = a(—y), b(y) = b(—y), a'(y) = —a!(—y) and 
b'(y) = ~b'(—y)), then, in the limit e — > 0, one gets 



a 
a 



\y=0 



\abp. 



Similarly, integrating the r/'-component and using the last equation gives 



b' 



\y=0 



--ab{p + p). 



(21) 



(22) 



These two conditions are called the junction conditions and play an important role in 
cosmology, describing how the brane with some given energy-momentum tensor 
can be embedded in a higher-dimensional space-time. The other components of the 
Einstein equations should be compatible with these conditions. It is not difficult to 
show that the restriction of the 05 component to y = leads to 



p + 3-(p + p) =0, 



(23) 



where we have made use of the junction conditions ()21j) and (|22|) . This is nothing 
but matter conservation on the brane, which must hold because we have assumed that 
matter is confined on the brane. 

Restricting the 55-component of the bulk equations to the brane and using again 
the junction conditions gives 



a ab , , , 

+ kb 2 

a ab 



« 2&2 r 1 , 



(24) 



Changing to cosmic time dr 
conservation gives ([SI], [3*5] ) 

d(H 2 e ia ) 2 



abdt, writing a = exp(a(t)) and using the energy 



A 5 e 4Q + 



d 



P 



da 3 da \ 36/ 

In this equation aH = da/ dr. This equation can easily be integrated to give 



H 2 



A 5 a 



36 6 



(25) 



(26) 



Let us split the total energy density and pressure into parts coming from matter and 
brane tension, i.e. we write p = pu + o~ and p = pu — o~. Then we obtain Friedmann's 
equation 



H 2 



8ttG 



-PM 



1 + 



Pm 
2a 



A 4 a 
+ tt + 4 



where we have made the identification 



(27) 

(28) 
(29) 

Imposing the fine-tuning relation ("T""]) of the static Randall-Sundrum solution in the last 
equation, we see that A 4 = 0. If there is a small mismatch between the brane tension and 



8nG 


a 




3 


" 18 




A 4 


_ a 2 


As 


3 


" 36 


- — 

6 
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the five-dimensional cosmo logical constant, an effective four-dimensional cosmo logical 
constant is generated. Hence, the Randall-Sundrum setup does not provide a solution 
to the cosmological constant problem, as one has to impose a relation between brane 
tension and the cosmological constant in the bulk, i.e. equation (|T4*|) . 

Another result is that the four-dimensional Newton constant is directly related to 
the brane tension in this model. 

The constant /i appears in the derivation above as an integration constant. The 
term including p is sometimes called the dark radiation term (see e.g. [SE]-|SH1)- The 
parameter p can be obtained from a full analysis of the bulk equations jSHl-jll] (we will 
discuss this in section 4). An extended version of Birkhoff's theorem tells us that, if 
the bulk space-time is AdS, this constant is zero |42| . If the bulk is AdS-Schwarzschild 
instead, p is non-zero but a measure of the mass of the bulk black hole. In the following 
we will assume that p = and A4 = 0. 

The most significant change in Friedmann's equation compared to the usual four- 
dimensional form is the appearance of a term proportional to p 2 . This term implies that 
if the matter energy density is much larger than the brane tension, i.e. pu ^> cr, the 
expansion rate H is proportional to pm, instead of ^Pm- The expansion rate is, in this 
regime, larger in this brane world scenario. Only in the limit where the brane tension is 
much larger than the matter energy density is the usual behaviour H oc ^fp~u recovered. 
This is the most important change in brane world scenarios. As we will see later, it is 
quite generic and not restricted to the Randall-Sundrum brane world model. 

At the time of nucleosythesis the brane world corrections in Friedmann's equation 
must be negligible, otherwise the expansion rate is modified and the results for the 
abundances of the light elements are unacceptably changed. This implies that a > 
(lMeV) 4 . Note, however, that a much stronger constraint arises from current tests for 
deviation from Newton's law [13] (assuming the Randall-Sundrum fine-tuning relation 
()14j) ): Kg 3 > 10 5 TeV and a > (lOOGeV) 4 . Similarily, cosmology constrains the amount 
of dark radiation. It has been shown that the energy density in dark radiation can be 
at most be 10 percent of the energy density in photons [3%] . 

Finally we can derive Raychaudhuri's equation from Friedmann's equation and the 
energy-conservation equation: 



2.1.2. Another derivation of Einstein's equation In the following, we present a more 
geometrical and, hence, more powerful derivation of Einstein's equation on the brane 



For this, consider an arbitrary (3+1) dimensional hypersurface M. with unit normal 
vector n a embedded in a five-dimensional space-time. The induced metric and the 
extrinsic curvature of the hypersurface are defined as 




dH 



(30) 



h\ =5 a b -n a n b , 
K ab = h a c h b d V c n d . 



(31) 
(32) 
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For the derivation we need three equations; two of them relate four-dimensional 
quantities constructed from h ab to full five-dimensional quantities constructed from 
g a b- We just state these equations here and refer to for their derivation. The first 
equation is the Gauss equation, which reads 

R { ah = hJh b k h c l h d m R jk i m - 2K a[c K d]b . (33) 

This equation relates the four-dimensional curvature tensor R a 4 bcd i constructed from h a b, 
to the five-dimensional one and K a b- The next equation is the Codazzi equation, which 
relates K ab , n a and the five-dimensional Ricci tensor: 

Vi 4) n - V®K = rfh\R hc . (34) 

One decomposes the five-dimensional curvature tensor R a bcd into the Weyl tensor C a bcd 
and the Ricci tensor: 

2 1 

Rabcd = g (g a [cRd}b — gb[cRd]a^j — -^Rg a [ c gb}d + C a bcd- (35) 

If one substitutes the last equation into the Gauss equation and constructs the four- 
dimensional Einstein tensor, one obtains 

= I (G cd h c a h d b + (G cd n c n d - Ig) h ab 



+ KK ah - K a c K bc -\(K 2 - K cd K cd ) h ab - E ab , (36) 

where 

E ab = C ac b d n c n d . (37) 

We would like to emphasise that this equation holds for any hypersurface. If one 
considers a hypersurface with energy momentum tensor T a b, then there exists a 
relationship between K a b and T a b {T is the trace of T a b) 



[K ab ] = -4 (Tab - -habTj , (38) 

where [...] denotes the jump: 

[f](y) = \im^ Q (f(y + e)-f(y-e)). (39) 

These equations are called junction conditions and are equivalent, in the cosmo logical 
context, to the junction conditions (|2T|) and Splitting T a b = T~ab—o~h a b and inserting 
the junction condition into equation (|36p. we obtain Einstein's equation on the brane: 

G$ = 8nGr ab - A 4 h ab + 4p ab - E ab . (40) 

The tensor 7i ab is defined as 

TTafe = -TT^Tab - -Jac^ + ■^h, oh T cd T cd - —T 2 h ab , (41) 

iZ 4 o Z4 



whereas 



K 



1 



SttG = '^a (42) 
6 



A 4 =^ A 5 + > 2 . (43) 
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Note that in the Randall-Sundrum case we have A 4 = due to the fine-tuning between 
the brane tension and the bulk cosmo logical constant. Moreover E ab = as the Weyl- 
tensor vanishes for an AdS space-time. In general, energy conservation and the Bianchi 
identities imply that 



on the brane. 

Clearly, this method is very powerful, as it does not assume homogeneity and 
isotropy nor does it assume the bulk to be AdS. In the case of an AdS bulk and a 
Friedmann-Robertson- Walker brane, the previous equations reduce to the Friedmann 
equation and Raychaudhuri equation derived earlier. 

To conclude this section, we have given two derivations of Einstein's equation as 
perceived by an observer confined on a brane embedded in an Anti-de Sitter bulk. In the 
derivation the emphasis was given to the brane observer. This point of view, however, 
neglects the role of the bulk. In particular it does not say anything about the magnitude 
of E a i) once the brane becomes inhomogeneous. In fact, the set of equations on the 
brane are not closed in general This is in particular a problem when discussing 
the evolution of cosmological perturbations: the evolution of the bulk gravitational field 
can not be solved using Einstein's equation on the brane alone, needing instead a full 
five-dimensional analysis. We will come back to this point in a later section. 

2.2. Slow-roll inflation on the brane 

We have seen that the Friedmann equation on a brane is drastically modified at high 
energy where the p 2 term dominates. As a result the early universe cosmology in brane 
world models tends to be different from that of standard four-dimensional cosmology. 
Thus, it seems natural to look for brane effects on early universe phenomena such as 
inflation (see in particular [IE] and [1H]) and on phase-transitions [HTT] . 

Let us therefore consider a scalar field confined on the brane. The energy density 
and the pressure of this field are given by 



where V(<f>) is the potential energy of the scalar field. The full evolution of the scalar 
field is described by the (modified) Friedmann equation, the Klein-Gordon equation 
and the Raychaudhuri equation. 

2.2.1. Single field slow-roll inflation In the following we will discuss the case of slow- 
roll inflation. Therefore, the evolution of the fields is governed by (from now on a dot 
stands for the derivative with respect to cosmic time) 



4V a Kab = V a E ab 



(44) 



(45) 



(46) 




(47) 
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It is easy to show that these equations imply that the slow-roll parameters are given by 

(49) 



* = 1 




2 


~4a(a 


+ V)] 


H 2 16nG 




) 


_(2a- 


-V)\ 



V" 

V = 



3H 2 8ttG 




2a 



(50) 



2a + V} 

The modifications to General Relativity are contained in the square brackets of these 
expressions. They imply that for a given potential and given initial conditions for the 
scalar field the slow-roll parameters are suppressed compared to the predictions made in 
General Relativity. In other words, brane world effects ease slow-roll inflation [48J. In 
the limit a <^V the parameters are heavily suppressed. It implies that steeper potentials 
can be used to drive slow-roll inflation [49J. What does this mean for cosmological 
perturbations generated during inflation? In order to study cosmological perturbations, 
one needs to have information about the projected Weyl tensor, which encodes the 
dynamics of the bulk space-time. As already said, one needs a full five-dimensional 
analysis to study the evolution of SE^ U . However, one can get some information by 
neglecting the backreaction of the bulk, i.e. by neglecting the contribution of E^. Let 
us then consider scalar perturbations first. The perturbed line element on the brane has 
the usual four-dimensional form 

ds 2 = -(1 + 2A)dt 2 + 2d i Bdtdx i + ((1 - 2^)5^ + D ij E)dx i dx j , (51) 

where the functions A, B, E and i/j depend on t and x % . We will make use of the gauge 
invariant quantity [22| 

C = ^ + H 5 -P. (52) 
P 

In General Relativity, the evolution equation for ( can be obtained from the energy 
conservation equation [22], which reads on large scales 



H 

C = ; — <Wd- (53) 

P + P 

In this equation <5p na d = Sp to t — c 2 s 8p is the non-adiabatic pressure perturbation. Thus, 
in the absence of non-adiabatic pressure perturbations the quantity ( is conserved on 
superhorizon scales. The energy conservation equation holds for the Randall-Sundrum 
model as well. Therefore, eq. ()53|) is still valid for the brane world model we consider. 
For inflation driven by a single scalar field, 5p na d vanishes and therefore ( is constant on 
superhorizon scales during inflation. Its amplitude is given in terms of the fluctuations 
in the scalar field on spatially flat hypersurfaces: 

The quantum fluctuation in the (slow-rolling) scalar field obeys ((5(f)) 2 ) ~ (H/2ti) 2 , as 
the Klein-Gordon equation is not modified in the brane world model we consider. The 
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amplitude of scalar perturbations is [HI] A 2 S = 4(£ 2 )/25. Using the slow-roll equations 
and eq. one obtains 



Al 




( 5127T 



\75M$ 



\k=aH 



(55) 



Again, the corrections are contained in the terms in the square brackets. For a given 
potential the amplitude of scalar perturbations is enhanced compared to the prediction 
of General Relativity. 

The arguments presented so far suggest that, at least for scalar perturbations, 
perturbations in the bulk space-time are not important during inflation. This, however, 
might not be true for tensor perturbations, as gravitational waves can propagate into 
the bulk. For tensor perturbations, a wave equation for a single variable can be derived 
[53] . The wave equation can be separated into a four-dimensional and a five-dimensonal 
part, so that the solution has the form h%j = A(y)h(x fl )eij, where is a (constant) 
polarisation tensor. One finds that the amplitude for the zero mode tensor perturbation 
is given by [EH] 

' -H 2 F 2 (H/n)\ k=aII , 



with 



F(x) 



\/l + x 2 — x 2 sinh 1 ( — 

x 



-1/2 



where we have defined 



1/2 



pi- 



(56) 



(57) 



(5f 



\4iraJ 

It can be shown that modes with m > 3H/2 are generated but they decay during 
inflation. Thus, in this scenario one expects only the massless modes to survive until the 
end of inflation [EE] , [EH] . From eqns. ([ST?]) and (jB3]l one sees that the amplitudes of scalar 
and tensor perturbations are enhanced at high energies, though scalar perturbations are 
more enhanced than tensors. Thus, the relative contribution of tensor perturbations 
will be suppressed if inflation is driven at high energies. 

Interestingly, to first order in the slow-roll parameters, the consistency relation 
between the amplitudes of scalar and tensor spectra and the spectral index of the tensor 
power spectrum 



(59) 



found in General Relativity holds also for inflation driven on the brane in AdS |57j . 
However, it was shown in |3HJ- that this degeneracy is broken at second order in 
the slow-roll parameters. Nevertheless, it will be difficult to distinguish observationally 
between models of inflation in General Relativity and the models of brane world inflation 
considered here. For a recent discussion on constraints on brane world inflation, see 
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2.2.2. Two field slow-roll inflation Let us consider now the case for two slowly rolling 
scalar fields, and \i driving a period of inflation in the brane world model we consider 
[B4*j . To ensure the slow roll behaviour of the fields, the slow-roll parameters 
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(60) 



should be small. However, instead of working with the fields and \ directly, it is 
more convenient to consider the so-called adiabatic a and entropy field s. Fluctuations 
in s are related to entropy perturbations, whereas perturbations in a describe purely 

The field a and s are defined through a field rotation: 



adiabatic perturbations 
5a = (cos 
6s = — (sin 

with 

cos 9 = 



+ (sin 9)5x 
+ (cos6)5x, 

■. , sin 9 



X 



X 2 



X 2 



The slow-roll parameters for o and s are given by 



16ttG 



Vo 

V 



77 + e 



x- 



(61) 
(62) 

(63) 
(64) 



and 



Vo-a = V<t>4> cos2 + ^V<t>x cos 9 sin 9 + r] xx sin 2 9, 

Vss = V&<f> srn2 ~~ 2?y^ x cos # sin + r/ xx cos 2 0, (65) 

^<7 S = (?7 X x ~ W) sin6»cos6» + r^ x (cos 2 6> - sin 2 6»). 

It can be shown that the total curvature perturbation ( and the entropy 
perturbation 

HSs 



S 



a 



(66) 



obey the equations 

C ~ - 2H Vas S (67) 

5~ [-2e-ris a + V^]HS. (68) 

These are the same expressions as found in General Relativity [HS]. Note, however, 
that at very high energies, i.e. V ^> A, we have that the ratios (/S and S/S are 
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suppressed by a factor « y\/V, when compared to General Relativity. It implies that, 
for given potentials and initial conditions, entropy perturbations have less influence 
on adiabatic perturbations. Even if entropy perturbations are generated, the final 
correlation between the adiabatic and entropy modes would be suppressed. 



2.3. Primordial black holes on the brane 

Let us now turn our attention to another interesting consequence of brane worlds: the 
modification of the evaporation law of primordial black holes, which results in a higher 
lifetime and lower temperature at evaporation jSS],[nZI- It was shown that the mass- 
temperature relation for a black hole confined on the brane with radius ro much smaller 
than the AdS curvature length scale is modified. Consider such a small black hole on 
the brane. Near the horizon, the metric is (approximately) that of a 5D Schwarzschild 
black hole, which reads, using spherical coordinates (/ = 1 — r^/r 2 ) 

dsj = -f(r)dt 2 + / _1 (r)dr 2 + r 2 dVt\. (69) 

The near-horizon geometry induced on the brane is given by 

ds\ = -f(r)dt 2 + f-\r)dr 2 + r 2 dQ 2 2 , (70) 

which is different from the four-dimensional Schwarzschild metric. However, as in 
General Relativity, it is the vicinity of the horizon which determines the thermodynamic 
quantities such as entropy and temperature. For the temperature one finds 



'Mr 



1/2 



t - = ^(ifI • < 71) 



Note that Tbh oc 1/vM instead of Tbh oc 1/M as in General Relativity. Furthermore, 
it can be shown that Stefan's law for a D-dimensional black hole reads jHZ| 

^ a g D o- D A e{i , D T D , (72) 
dt 

in which gjj denotes the total number of effective bosonic and fermionic degrees of 
freedom, o~d is the D-dimensional Stefan-Boltzmann constant defined per degree of 
freedom (C{D) is the Riemann zeta function and Qd-2 the volume of a D — 2 sphere) 

"> = ^hnD)m (73) 

and A e g t -£) is the effective black hole area calculated from the induced metric. 

This implies that, for a given mass, the lifetime t cvap ~ M/M of such a small black 
hole on the brane can be much longer than black holes in General Relativity. One finds 
that, for a small five-dimensional black hole, 

(i\ 2 

^evap,5D ~ ^evap,4D — , (74) 
W 

where t ev ap,4D is the lifetime of a black hole of the same mass as predicted with General 
Relativity. In fact, in 67 ,[68j it was shown that such a small black hole can survive 
until the present epoch. 
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Once a small black hole has formed, it will accrete matter and radiation surrounding 
it. It is thought that in the standard cosmology this process does not lead to a significant 
enhancement of the black hole mass [HE]- However, this is no longer true in brane world 
cosmology in the high energy regime, where Friedmann's equation is modified. In this 
regime the radiation density obeys p ra( j oc t -1 instead of p ra d oc t~ 2 . As a result, the 
growth in the mass of the black hole is M oc t^ 1 (in General Relativity one finds 
M oc t~ 2 ). As discussed in depth in [HS], the consequence of this is that a small black 
hole can experience significant growth in the high energy regime, albeit with details 
depending on the accretion efficiency Once the high energy regime is over, the accretion 
history is the same as in standard cosmology. In 7Q 1 one finds a detailed discussion on 
how these considerations modify some of the known constraints on primordial black 
holes. 

For more work on black holes in brane cosmology, see [2]- [HI]. For reviews, see 
[HI and HU. 



2-4- Graviton Backreaction 

We now discuss the production of gravitons by processes on the brane [23, [lEj- As 
branes carry matter fields and gravity couples to matter, one can envisage the radiation 
of gravitons into the bulk. The effect of the graviton emission is to modify the 
background geometry. In turn, this leads to the expected dark radiation term on the 
brane. As a first step we will neglect the back-reaction of gravitons impinging on the 
brane. This is much more involved and will be briefly touched upon at the end of the 
subsection. 



2.4-1- Dark Radiation Physically, matter particles on the brane can collide and lead to 
the presence of bulk gravitons. As such the brane gravitons carry an energy-momentum 
tensor and therefore lead to a modification of the bulk Einstein equations. Let us first 
assume that the gravitons escape radially from the brane. The energy-momentum tensor 
in the bulk is then 

Tab = o-k a k b (75) 

where k a is the null impulsion of the gravitons. The general situation when radiation is 
not radial will be dealt with later. One can choose the ansatz 

ds 2 = -f(r, v)dv 2 + 2drdv + r 2 dx 2 (76) 

where 

f(r,v) = k 2 r 2 -^- (77) 

for the bulk metric. This is the so-called Vaidya metric. The function C(v) generalises 
the black hole mass as it occurs in the static solutions obeying Birkhoff 's theorem. The 
Weyl parameter C(v) satisfies 

dC 2k 2 a 3 /. / , ; , 9 \ 2 

_ = _p> < T - 4J^7A (78) 
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in terms of the proper time on the brane. The brane junction conditions are now 
modified and lead to (H — f/r) 



This is the Friedmann equation where the Weyl parameter C(v) plays the role of the 
dark radiation term. Similarly the total matter and pressure contents, p m = A + p and 
Vm = —A + p, satisfy a non-conservation equation 



r 



p m + 3- (p m + p m ) = -2a (80) 
r 

One can see that the effect of the graviton emission is to induce a loss term in the 
conservation equation. The cosmology on the brane is now highly modified by the 
graviton radiation term. To go further we need to evaluate a. 

Gravitons in the bulk are produced according to the tree level process ip + ^ — > 
graviton where ip and ip are particles and antiparticles confined to the brane. As long 
as the temperature on the brane is much lower than both re^ 3 and re J 1 , then the 
temperature T is always bigger than the Hubble rate. Indeed at high energy H ~ re 2 T 4 
implies that H T and similarly at low energy H ~ k^T 2 implying that H T 
in this regime as well. We have used Stefan's law whereby p ~ T 4 in the radiation 
epoch. Hence cosmological effects can be neglected and the scattering cross sections can 
be computed in Minkowski space on the brane. The interaction Lagrangian between 
matter fields and gravitons is given by 

S = k 5 J dmu m (0) J d 4 XT^hffi (81) 

where m labels the continuous spectrum of Kaluza-Klein states. The wave function of 
the gravitons in the extra dimension is u m (y) corresponding to the four dimensional field 
hffi. Here r^ v is the energy momentum of particles on the brane. The matrix element 
corresponding to this tree level interaction can be computed (averaged over spins) 



2 | • • b 



\M\ = 4\u m {Q)\A- (82) 

where s is the Mandelstam variable s = (pi + P2) 2 corresponding to the annihilation of 
particles with momenta pi j2 . The factor A is given by A = 2/3 for scalars, A = 1 for 
fermions and A = 4 for photons. 

The Boltzmann equation in curved space is 

p + 3H(p + p) = ~JdmJ T0^C m (83) 

where the collision term reads 

^ = y§y^\^f lf ^( Pm - Pl - P2) (84) 

where /j = l/(e E ^ T ± 1) depending on the statitics of the particles. Here p m is the 
momentum of the gravitons. Neglecting the contribution of the light gravitons m k, 
and using u m (0) = l/y^ for m ^> k one finds that 

' 315C(9/2)C(7/2) ~ 

5127T 2 



p + AHp = - 



g(T)4T* (85) 
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where the effective number of species is given by 

g(T) = ^ + Ag v + (1 - 2- 7 / 2 )(l - 2^)3 f (86) 

Here g s , gf and g v are the number of relativistic particles of spin 0, 1/2 and 1. 
This allows one to identify the parameter a 

o = ^4P 2 (87) 

where Stefan's law reads p = ^^T 4 where g* = g s + 9v + 7/&9f an d we have defined 
the parameter 

212625 w . xw . .g(T) , s 

" = -64^ C(9/2)C(7/2) ^T (88) 
The cosmological equations can now be solved in the high energy regime leading to the 
Weyl parameter on the brane 

C(t) ~ aa 4 (89) 

and 

a(t) ~ t 1 /^ (90) 

Numerically a ~ 0.02 in the standard model of particle physics. The correction to the 
t 1//4 behaviour in the high energy phase is tiny. Nevertheless, notice that the ratio 

ew = >~r (91) 

Prad 4 

at low energy. The presence of some extra radiation energy density during 
nucleosynthesis is highly constrained leading to 

e w < 8.10~ 2 (92) 
which is marginally higher than a/4. Forthcoming constraints from CMB anisotropics 
will lead to a tightening of this bound. This might open up the possibility of confronting 
the prediction of brane models with experiment. 

In the following we will present the case when radiation is not radial. 



2.4-2. Gravitons bouncing off the brane So far it has been assumed that radiation is 
radial. This is a drastic assumption. In particular this implies that gravitons cannot 
be emitted and then captured by the brane. An analysis which includes the effect of 
gravitons bouncing off the brane has been performed in [ZB]- The gravitons emitted 
after a particle-antiparticle annihilation on the brane are not all radial. In particular, 
the energy-momentum tensor is not well approximated by the form given in (|75J) but 
also contains non-radial terms proportional to p 2 . Moreover, the gravitons may follow 
geodesies in the bulk which could lead to subsequent collisions with the brane. As a 
result, the non-conservation equation for the dark component depends on the amount of 
radiation being emitted and bouncing off the brane. This has been studied numerically. 
In particular the ratio ew has been evaluated and is enhanced in comparison with the 
value, a/4, given in the previous section. Consequently, this tightens the bounds coming 
from nucleosynthesis and CMB. Finally, the dynamical consequences of a non-constant 
Weyl parameter in radiating branes have been investigated in [77] . 
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2.5. Non-Z 2 Symmetric Branes 

Most brane world scenarios assume a Z 2 symmetry about our brane, motivated by M- 
theory. However, a lot of recent models are not directly derived from M-theory and 
the motivation for maintaining the Z 2 symmetry is less clear. For example, multi-brane 
scenarios have been suggested which involve branes that, although lying between two 
bulk space times with the same cosmological constant, do not possess a Z 2 symmetry 
of the metric itself [78\. In this subsection we investigate the cosmological implications 
of relaxing the Z 2 symmetry about the brane world embedded in AdS-Schwarzschild 

mm- 

When one does not impose a Z 2 symmetry about our brane world the Friedmann 
equation on the brane acquires an extra term for the lack of Z 2 symmetry The method 
of deriving the Friedmann equation is as outlined previously. However, the junction 
conditions are modified. Let us write [HO] 

a'(p + ) = -a'(p-) + d(t). (93) 

where d(t) is some function of time only and has yet to be determined. d(t) represents 
the asymmetry of the metric across the brane and the fact that it is not identically zero 
for all time is what is referred to as the brane being i non-Z 2 symmetric' in the rest of 
this review. Calculating the Israel junction conditions one finds that 

2F 

d(t) = —j, (94) 
p b a 6 

where F is an integration constant which, when non-zero, dictates to what extent the 
Z 2 symmetry is broken (in fact, it is related to the difference of the black hole masses 
on each sides of the branes, see [22] )■ This results in a modified Friedmann equation 

d\ 2 k 4 pa k 4 2 k C F 2 , . 



V 18 r 36 r o? a 4 (p + p x ) 2 a 8 ' 
where p\ is the energy density associated to the brane tension. So the absence of the 
Z 2 symmetry gives rise to an extra term in the Friedmann equation. For a radiation 
dominated Universe where p = 7/a 4 , the extra term behaves as F 2 / r y 2 as p — > 00. 
This results in the Hubble constant being approximately constant at very early times, 
resulting in an inflation-like period. At late times the extra term behaves as (Fp/ / yp\) 2 
as p — > 0, so the effects are no longer significant and our solution reverts to that analysed 
previously. 

The effect of this extra term in the Friedmann equation has been analysed in |8()j . 
where it was shown that the universe undergoes a period of exponential expansion during 
a radiation dominated epoch. However, given that the universe must be described by 
the usual Friedmann equation by the time of nucleosynthesis the maximum expansion 
due to the F-term was found to be 10 4 , which is less than that needed to solve the 
standard cosmological problems. 
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2. 6. Some final remarks 

The Randall-Sundrum brane world is the simplest brane world model, which leads to a 
consistent cosmological framework. However, the model is related to deeper issues, such 
as the AdS-CFT correspondence (see e.g. [HT]-|86j) and therefore to the holographic 
principle [HZ|. We cannot go into the details here. Similarly, we do not cover inflation 
driven by the trace anomaly of the conformal field theory living on the brane (see e.g. 
[HH]-|ni])- These are interesting and important developments which give insights into 
gravity and the consequences for early universe cosmology. 

3. Including a Bulk Scalar Field 

In this section we are going to generalise the results found in the previous section where 
the bulk was empty to a more realistic case where the effects of a bulk scalar field 
are included. This is motivated by supersymmetric extensions to the Randall-Sundrum 
model, and of course is a more realistic limit of string theory which automatically has the 
string dilaton field. The projective approach considered previously, where one focuses 
on the dynamics of the brane, can be extended to this case. Here one studies both 
the projected Einstein and the Klein-Gordon equation [OB]- Again, bulk effects 
do not decouple, so the dynamics are not closed. As we will see, there are now two 
objects representing the bulk back-reaction: the projected Weyl tensor E^ u and the loss 
parameter A$ 2 . In the case of homogeneous and isotropic cosmology on the brane, the 
projected Weyl tensor is determined entirely up to a dark radiation term. Unfortunately, 
no information on the loss parameter is available. This prevents a rigorous treatment 
of brane cosmology in the projective approach. However, the inclusion of a bulk scalar 
field opens up other avenues. As we will see later in this section, the bulk scalar field 
can give rise to late-time acceleration of the universe, i.e. it could act as a quintessence 
field, and also leads to a variation in Newton's constant. 

Another route to analysing brane models is to study the motion of a brane in a 
bulk space-time. This approach is successful in the Randall-Sundrum case thanks to 
Birkhoff's theorem which dictates a unique form for the metric in the bulk ^2] . In the 
case of a bulk scalar field, no such theorem is available. One has to resort to various 
ansatze for particular classes of bulk and brane scalar potentials (see e.g. 
We discuss this in more detail in section 4. 

3. 1 . BPS Backgrounds 

3.1.1. Properties of BPS Backgrounds As the physics of branes with bulk scalar fields 
is pretty complicated, we will start with a particular example where both the bulk and 
the brane dynamics are related and under control (see also |l()4j and |l()5j ). This is 
the situation in the BPS case and is automatic in the supersymmetric case. We specify 
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the bulk Lagrangian as 

S = / <*V=SS (R - \ {{d<pf + V{4>))) (96) 

where V(<p) is the bulk potential. The boundary action depends on a brane potential 

U B (<f>) 

3 J g?W=^b(0o) (97) 



Sb - - — 
2k 5 



where Ub(4>o) is evaluated on the brane. The BPS backgrounds arise as a special case of 
this general setting with a particular relationship between the bulk and brane potentials. 
This relation appears in the study of N = 2 supergravity with vector multiplets in the 
bulk. The bulk potential is given by 



'dW\ 



2 



V = [ —J - W 2 (9? 



where W(<f>) is the superpotential. The brane potential is simply given by the 
superpotential 

U B = W (99) 

We would like to mention that the last two relations have also been used to generate 
bulk solutions without necessarily imposing supersymmetry [9~o]. |lU2| . Notice that the 
Randall-Sundrum case can be retrieved by putting W = est. Supergravity puts further 
constraints on the superpotential which turns out to be of the exponential type |lU3j 

W = Ake^ (100) 

with a = — 1/ vT2, 1 / v3- I n the following we will often choose this exponential potential 
with an arbitrary a as an example. The actual value of a does not play any role and 
will be considered generic. 

The bulk equations of motion comprise the Einstein equations and the Klein-Gordon 
equation. In the BPS case and using the following ansatz for the metric 

ds 2 = a{yfr] lxv dx^dx 1 ' + dy 2 , (101) 

these second order differential equations reduce to a system of two first-order differential 
equations 

Notice that when W = est one easily retrieves the exponential profile of the Randall- 
Sundrum model. 

An interesting property of BPS systems can be deduced from the study of the 
boundary conditions. The Israel junction conditions reduce to 
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and for the scalar field 

4>\B = -7Tt\ 



(104) 



This is the main property of BPS systems: the boundary conditions coincide with 
the bulk equations, i.e. as soon as the bulk equations are solved one can locate the 
BPS branes anywhere in this background, there is no obstruction due to the boundary 
conditions. In particular two-brane systems with two boundary BPS branes admit 
moduli corresponding to massless deformations of the background. They are identified 
with the positions of the branes in the BPS background. We will come back to this later 
in section 5. 

Let us treat the example of the exponential superpotential. The solution for the 
scale factor reads 

a = (1 - Aka 2 y) 1/4a \ (105) 
and the scalar field is given by 

= — Wl -Aka 2 y). (106) 
a 

For a — > 0, the bulk scalar field decouples and these expressions reduce to the Randall- 
Sundrum case. Notice a new feature here, namely the existence of singularities in the 
bulk, corresponding to 

a(v)l* = (107) 
In order to analyse singularities it is convenient to use conformal coordinates 

du = (108) 

In these coordinates light follows straight lines u = ±t. It is easy to see that the 
singularities fall into two categories depending on a. For a 2 < 1/4 the singularity 
is at infinity = oo. This singularity is null and absorbs incoming gravitons. For 
a 2 > 1/4 the singularity is at finite distance. It is time-like and not wave- regular, i.e. 
the propagation of wave packets is not uniquely defined in the vicinity of the singularity. 
For all these reasons these naked singularities in the bulk are a major drawback of brane 
models with bulk scalar fields 106 . In the two-brane case the second brane has to sit in 
front of the naked singularity. This is automatic in the supersymmetric required 
by conservation of global charges. 

3.1.2. de Sitter and anti de Sitter Branes Let us modify slightly the BPS setting by 
detuning the tension of the BPS brane. This corresponds to adding or subtracting some 
tension compared to the BPS case 

U B = TW (109) 

where T is a real number. Notice that this modification only affects the boundary 
conditions; the bulk geometry and scalar field are still solutions of the BPS equations 
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of motion. In this sort of situation, one can show that the brane does not stay static. 
For the detuned case, the result is either a boosted brane or a rotated brane. Defining 
by u(x^) the position of the brane in conformal coordinates, one obtains 

{duf = (110) 

The brane velocity vector d^u is of constant norm. For T > 1, the brane velocity vector 
is time-like and the brane moves at constant speed. For T < 1 the brane velocity vector 
is space-like and the brane is rotated fl07[ HU8j . Going back to a static brane, we see 
that the bulk geometry and scalar field become dependent. In the next section we 
will find many more cases where branes move in a static bulk or equivalently, a static 
brane borders a non-static bulk. 

Let us now conclude this section by studying the brane geometry when T > 1. In 
particular one can study the Friedmann equation for the induced scale factor 

H 2 = 7 ^W\ (111) 
lo 

where W is evaluated on the brane. Of course we find that cosmological solutions are 
only valid for T > 1. Now in the Randall- Sundrum case W = 4k leading to 

H 2 = (T 2 - l)k 2 . (112) 

In the case T > 1 the brane geometry is driven by a positive cosmological constant. This 
is a de Sitter brane. When T < 1 the cosmological constant is negative, corresponding 
to an AdS brane. We are going to generalise these results by considering the projective 
approach to the brane dynamics. 

3.2. Bulk Scalar Fields and the Projective Approach 

3.2.1. The Friedmann Equation We will first follow the traditional coordinate- 
dependent path. This will allow us to derive the matter conservation equation, the 
Klein-Gordon and the Friedmann equations on the brane. Then we will concentrate on 
the more geometric formulation where the role of the projected Weyl tensor will become 
transparent |114| . [TTK] . Again, in this subsection we will put k§ = 1. 

We consider a static brane that we choose to put at the origin y = 0. and impose 
6(0, t) = 1. This guarantees that the brane and bulk expansion rates 

AH = d T ^\o, 3H B = d rV /=^| (113) 

coincide. We have identified the brane cosmic time dr = ab\odt. We will denote by 
prime the normal derivative d n = -^\od y . Moreover we now allow for some matter to be 
present on the brane 

7? ^ = (-Pm,Pm,P m ,P m )- (H4) 

The bulk energy-momentum tensor reads 

T ab = - (da^) - \g ab ((«90) 2 + V) . (115) 
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The total matter density and pressure on the brane are given by 

3 3 

P = Pm + ~U B , P = Pm- ~U B - (116) 

The boundary condition for the scalar field is unchanged by the presence of matter on 
the brane. 

The (05) Einstein equation leads to matter conservation 

p m = -3H(p m + p m ). (117) 
By restricting the (55) component of the Einstein equations we obtain 



H 2 = ^-^Q-^E + ^- (118) 



36 3 V 9 a 4 



in units of k\. The last term is the dark radiation, whose origin is similar to the 
Randall-Sundrum case. The quantities Q and E satisfy the differential equations 

Q + AHQ = HTl 
E + AHE = -pT 5 °. 

These equations can be easily integrated to give 



H 2 = & + ^ - — A I dA^ -2U)-—l draV^, (119) 
36 12 16a 4 J dr ' 12a 4 J h dr ' v ; 

up to a dark radiation term and we have identified 

^HtH- <i2o) 

This is the Friedmann equation for a brane coupled to a bulk scalar field. Notice that 
retarded effects springing from the whole history of the brane and scalar field dynamics 
are present. In the following section we will see that these retarded effects come from the 
projected Weyl tensor. They arise from the exchange between the brane and the bulk. 
Notice that Newton's constant depends on the value of the bulk scalar field evaluated 
on the brane (0o = <f>(t, y = 0)): 

8ttGW(0o) 4Ub(^o) 



(121) 

3 12 v ; 

On cosmological scales, time variation of the scalar field induces a time variation of 
Newton's constant. This is highly constrained experimentally |116j . [TT7] . leading to 
tight restrictions on the time dependence of the scalar field. 

To get a feeling of the physics involved in the Friedmann equation, it is convenient 
to assume that the scalar field is evolving slowly on the scale of the variation of the scale 
factor. Neglecting the evolution of Newton's constant, the Friedmann equation reduces 
to 

H P-+8-16 (122) 

Several comments are in order. First of all we have neglected the contribution due to 
the p 2 m term as we are considering energy scales below the brane tension. Then the 
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main effect of the scalar field dynamics is to involve the potential energy U and the 
kinetic energy 2 . Although the potential energy appears with a positive sign we find 
that the kinetic energy has a negative sign. For an observer on the brane this looks like 
a violation of unitarity. We will reanalyse this issue in section 5, when considering the 
low energy effective action in four dimensions and we will see that there is no unitarity 
problem in this theory. The minus sign for the kinetic energy is due to the fact that one 
does not work in the Einstein frame where Newton's constant does not vary, a similar 
minus sign appears also in the effective four-dimensional theory when working in the 
brane frame. 

The time dependence of the scalar field is determined by the Klein-Gordon equation. 
The dynamics is completely specified by 

4> + 4H<j> + -(- - Um)Pm-g^ = -qT + A$ 2 , (123) 
where p m = uj m p m . We have identified 

This cannot be set to zero and requires the knowledge of the scalar field in the vicinity 
of the brane. When we discuss cosmological solutions below, we will assume that this 
term is negligible. 

The evolution of the scalar field is driven by two effects. First of all, the scalar field 
couples to the trace of the energy-momentum tensor via the gradient of Ub- Secondly, 
the field is driven by the gradient of the potential U, which might not necessarily vanish. 

3.2.2. The Friedmann equation vs the projected Weyl tensor We are now coming back 
to the origin of the non-trivial Friedmann equation. Using the Gauss-Codazzi equation 
one can obtain the Einstein equation on the brane |92|.p3] 

Gab = -\uh ab + ^-r ab + n ab + \d a (t)d b (t) - -^-(d(j)) 2 h ab - E ab . (125) 
o 4 z lb 

Now the projected Weyl tensor can be determined in the homogeneous and isotropic 
cosmology case. Indeed only the Eqq component is independent. Using the Bianchi 
identity D a G a b = where D a is the brane covariant derivative, one obtains that 

' 12 , 3 TT \ 3 22 , 

leading to 



£ 00 + 4#£ o = d T + -Uj + -H^ + -f Pm (126) 



Eoo = ^J dra* U + |^ + + ^f Pm \ (127) 



Upon using 

Goo = 3# 2 (128) 

one obtains the Friedmann equation. It is remarkable that the retarded effects in the 
Friedmann equation all spring from the projected Weyl tensor. Hence the projected 
Weyl tensor proves to be much richer in the case of a bulk scalar field than in the empty 
bulk case. 
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3.2.3. Self-Tuning and Accelerated Cosmology The dynamics of the brane is not closed, 
it is an open system continuously exchanging energy with the bulk. This exchange is 
characterised by the dark radiation term and the loss parameter. Both require a detailed 
knowledge of the bulk dynamics. This is of course beyond the projective approach where 
only quantities on the brane are evaluated. In the following we will assume that the 
dark radiation term is absent and that the loss parameter is negligible. Furthermore, 
we will be interested in the effects of a bulk scalar field for late-time cosmology (i.e. 
well after nucleosynthesis) and not in the case for inflation driven by a bulk scalar field 



Let us consider the self-tuned scenario as a solution to the cosmological constant 
problem. It corresponds to the BPS superpotential with a — 1. In that case the potential 
U = for any value of the brane tension. The potential U — can be interpreted 
as a vanishing of the brane cosmological constant. The physical interpretation of 
the vanishing of the cosmological constant is that the brane tension curves the fifth 
dimensional space-time leaving a flat brane intact. Unfortunately, the description of the 
bulk geometry in that case has shown that there was a bulk singularity which needs to 
be hidden by a second brane whose tension is fine tuned with the first brane tension. 
This reintroduces a fine tuning in the putative solution to the cosmological constant 
problem |22j . 

Let us generalise the self-tuned case to a ^ 1, i.e. Ub = TW, T > 1 and W is the 
exponential superpotential. The resulting induced metric on the brane is of the FRW 
type with a scale factor 



leading to an acceleration parameter 



For the supergravity value a = — this leads to q = —4/7 and equation of state 
lo = —5/7. This is in coincidental agreement with the supernovae results. This model 
can serve as a brane quintessence model [103j. [TT4"] . We will comment on the drawbacks 
of this model later. See also |119j and |12Uj for similar ideas. 

3.2.4- The brane cosmological eras Let us now consider the possible cosmological 
scenarios with a bulk scalar field |114j . [Tl3] . We assume that the potential energy of 
the scalar field U is negligible throughout the radiation and matter eras before serving 
as quintessence in the recent past. 

At very high energy, i.e. energies above the tension of the brane, the non- 
conventional cosmology driven by the (? m term in the Friedmann equation is obtained. 
Assuming radiation domination, the scale factor behaves like 



(see e.g. |TU$j-[TLlj). 




(129) 




(131) 
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+ (3ln[^\ (132) 



and the scalar field 

^ci- 
In the radiation-dominated era no modification is present, provided 

= 0i (133) 

which is a solution of the Klein-Gordon equation as the trace of the energy-momentum 
tensor for radiation vanishes (together with a decaying solution, which we have 
neglected). In the matter-dominated era the scalar field evolves due to the coupling 
to the trace of the energy-momentum tensor. This has two consequences. Firstly, 
the kinetic energy of the scalar field starts contributing in the Friedmann equation. 
Secondly, the effective Newton constant does not remain constant. The cosmo logical 
evolution of Newton's constant is severely constrained since nucleosynthesis |llfi| . |ll7| . 
This restricts the possible time variation of <j). 

In order to be more quantitative let us come back to the exponential superpotential 
case with a detuning parameter T. The time dependence of the scalar field and scale 
factor become 

ft \§-ih 2 

a = ae [tJ 

where t e and a e are the time and scale factor at matter-radiation equality. Notice the 
slight discrepancy of the scale factor exponent with the standard model value of 2/3. 
The redshift dependence of the Newton constant is 

G N (z) = / z + l N 4 " 2 / 5 

G N ( Ze ) ~ Ue + J 

For the supergravity model with a = — and z e ~ 10 3 this leads to a decrease by 
(roughly) 37% since nucleosynthesis. This is marginally compatible with experiments 

Finally let us analyse the possibility of using the brane potential energy of the scalar 
field U as the source of acceleration now. We have seen that when matter is negligible 
on the brane, one can build brane quintessence models. We now require that this occurs 
only in the recent past. As can be expected, this leads to a fine-tuning problem as 

M 4 ~ Pc (135) 

where M 4 = (T — 1)§^ is the amount of detuned tension on the brane. Of course this is 

5 

nothing but a reformulation of the usual cosmological constant problem. Provided one 
accepts this fine tuning, as in most quintessence models, the exponential model with 
a = — is a cosmologically consistent quintessence model with a five-dimensional 
origin. 



(134) 



Brane World Cosmology 



28 



3. 3. Brief summary 

In this section we have explored the cosmology of a brane world model with a bulk 
scalar field, as motivated by supersymmetry. This is a complicated model, with more 
consequences than the Randall-Sundrum model. The main difference is that the 
gravitational constant becomes time dependent. As such it has much in common with 
scalar-tensor theories |118j . but there are important differences due to the projected 
Weyl tensor and its time evolution. The bulk scalar field can play the role of the 
quintessence field, as discussed above, but it could also play a role in an inflationary 
era in the very early universe (see e.g. |1U9| - [TT5] ). In any case, the cosmology of such 
a system is much richer and, because of the variation of the gravitational constant, 
more constrained. Similarly, it can leave a trace in the CMB anisotropies, which will 
also constrain the parameters of the theory, and result in the variation of fundamental 
constants. These results are discussed in sections 5 and 6 of this review. 

4. Birkhoff Theorem and its Violations 

We have already mentioned the extended Birkhoff theorem in section 2. It states that 
for the case of a vacuum bulk space-time, the bulk is necessarily static, in certain 
coordinates. A cosmologically evolving brane is then moving in that space-time, whereas 
for an observer confined to the brane the motion of the brane will be seen as an expanding 
(or contracting) universe. In the case of a scalar field in the bulk, a similar theorem 
is unfortunately not available, which makes the study of such systems much more 
complicated. We will now discuss these issues in some detail, following in particular 
jnn] and |lUUj . As we will see later in this section, the violation of Birkhoff theorem in 
the presence of a bulk scalar field is due to the projected Weyl tensor. 

4-1. Motion in AdS-Schwarzs child Bulk 

We have already discussed the static background associated with BPS configurations 
(including the Randall-Sundrum case) in the last section. Here we will focus on other 
backgrounds for which one can integrate the bulk equations of motion. Let us write the 
following ansatz for the metric 



where dT? is the metric on the three-dimensional symmetric space of curvature q = 
0, ±1. In general, the function A, B and R depend on the type of scalar field potential. 
This is to be contrasted with the case of a negative bulk cosmological constant where 
Birkhoff 's theorem states that the most general solution of the (bulk) Einstein equations 
is given by A 2 = f, B 2 = 1/ f and R = r where 



ds 2 = -A 2 (r)dt 2 + B 2 (r)dr 2 + R 2 {r)dT? 



(136) 



f{r) =q + 




(137) 
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We have denoted by I — 1/k — ^ —6 / (A^k 2 ) the AdS scale and [i the black hole mass 
(see section 2). This solution is the so-called AdS-Schwarzschild solution. 

Let us now study the motion of a brane of tension T/Hn such a background. The 
equation of motion is determined by the junction conditions. The resulting equation of 
motion for a boundary brane with a Z 2 symmetry is 

(r 2 + f(r)) 1/2 = jr (138) 

for a brane located at r j3Q|- Here r is the velocity of the brane measured with respect 
to the proper time on the brane. This leads to the following Friedmann equation 

^®' = ^-£ + $- im 

So the brane tension leads to an effective cosmological constant (T 2 — l)// 2 . The 
curvature gives the usual term familiar from standard cosmology while the last term 
is the dark radiation term whose origin springs from the presence of a black-hole in the 
bulk. At late time the dark radiation term is negligible for an expanding universe; we 
retrieve the cosmology of a FRW universe with a no n- vanishing cosmological constant. 
The case T = 1 corresponds of course to the Randall-Sundrum case. 

Notice that all the bulk physics in the absence of a bulk scalar field is captured by 
the unique choice of the static background solution. This solution is only parametrised 
by the AdS curvature and the black hole mass. These two numbers are all that is needed 
to describe the motion of any brane in a static bulk. We will now contrast this situation 
with the case of a bulk scalar field. 

4-2. Violating Birkhoff Theorem 

We now turn to a general analysis of the brane motion in a static bulk with a bulk scalar 
field. To do that it is convenient to parametrise the bulk metric slightly differently 

dr 2 

ds 2 = -f 2 (r)h(r)dt 2 + ttt + r 2 dY? a . (140) 

h[r) q 

where q = 0, ±1 is the spatial curvature. Now, the Einstein equations lead to (redefining 

*(*£ + v) (141) 

; , n 'A ( 142 ) 

and the Klein Gordon equation 

h( t>" + [y + Y + h '^' = ^- (143) 

Subtracting eq. (|141j) from (|142|) and solving the resulting differential equation, we 
obtain 

"drrcj)' 2 ]. (144) 
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This gives the link between the bulk scalar field and the bulk metric. 

Another piece of information is obtained by evaluating the spatial trace of the 
projected Weyl tensor. This is obtained by computing both the bulk Weyl tensor and 
the vector normal to the moving brane 

£. S -S = jif^y+x. (MM 

This is the analogue of the dark radiation term for a general background. The equations 
of motion can be cast in the form 

2 ley ^ 

f' ="f(/<-^-W 2 , (146) 



«' + ^-T (M ->° (147) 



K \V =Qn+-rH'-3^, (148) 



where we have defined 

H = (149) 

This allows us to retrieve easily some of the previous solutions. 
Choosing to be constant leads to / = 1, fi is constant and 

p r 

This is the AdS-Schwarzschild solution as specified by Birkhoff theorem. 
Another interesting case is obtained for spatially flat sections q = 

*V =~^T, (151) 



^ = "7^ + 4 (150) 



3 

d(- 

_2 



-d (r~ 4 ) , (152) 



K i v --^M) +H (l53) 

In this form it is easy to see that the dynamics of the bulk are completely integrable. 

First of all the solutions depend on an arbitrary function which determines the 
dynamics. Notice that 

/ = ^ (154) 

where no is an arbitrary constant. The radial coordinate r is obtained by simple 
integration of eq. ()151|) 

r = r e"^^ # . (155) 
Finally the rest of the metric follows from 

4k 2 , AM 



\V U?e 4 ^^ (156) 
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The potential V then follows from (|153|) . This is remarkable and shows why Birkhoff's 
theorem is not valid in the presence of a bulk scalar field. It is due to the infinite number 
of choices for the dark radiation function /i(0). Any choice of this function leads to a 
new static background in the bulk. 

Let us now discuss the cosmological consequences of this theorem. The Friedmann 
equation for a cosmological brane in the static background specified by fi{4>) is 

H 2 = n+^ 2 p 2 (157) 
3o 

where H is the Hubble parameter on the brane in cosmic time. In general the dynamics 
are extremely intricate. One can simplify and retrieve standard cosmology by studying 
the vicinity of a critical point ^# = 0. Parameterising 

6A 

fi = — + B<f) 2 (158) 



K 



5 

leads to the Friedmann equation 

# 2 = !(p 2 -^ 2 + 4+°(«- 4 ) ( l59 ) 

3d a 4 

Here 9 is an arbitrary integration constant. Notice that this is a small deviation from 
the Randall- Sundrum case as the scalar field behaves as 

= r~ B ' A (160) 

and goes to zero at large distances. Hence, standard cosmology is retrieved at low energy 
and long distance. Indeed, fi(4>) becomes constant and putting 6 = one obtains the 
brane Friedmann equation with its characteristic p 2 term. 



5. Cosmology of a Two— Brane System 

In this section we will once more include an ingredient suggested by particle physics 
theories, in particular M-theory. So far we have assumed that there is only one brane 
in the whole space-time. According to string theory, there should be at least another 
brane in the bulk. Indeed, in heterotic M-theory these branes are the boundaries of 
the bulk space-time [3]. Another motivation is the hierarchy problem. Randall and 
Sundrum proposed a two-brane model (one with positive and one with negative tension), 
embedded in a five-dimensional AdS space-time. In their scenario the standard model 
particles would be confined on the negative tension brane. As they have shown, in 
this case gravity is weak due to the warping of the bulk space-time. However, as will 
become clear from the results in this section, in order for this model to be consistent 
with gravitational experiments, the interbrane distance has to be fixed |122j . This can be 
achieved, for example, with a bulk scalar field. As shown in |122| and |123j . gravity in the 
two-brane model of Randall- Sundrum is described by a scalar-tensor theory, in which 
the interbrane distance, the radion, plays the role of a scalar field. Introduction of a bulk 
scalar field modifies the Brans-Dicke parameter (see |124j and |125j ) and introduces a 
second scalar field in the low-energy effective theory. Hence, in the case of two branes 
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and a bulk scalar field, the resulting theory at low energy is a bi-scalar-tensor theory 



In the following we will investigate the cosmological consequences when the distance 
between the branes is not fixed (for some aspects not covered here see e.g. |128| - [T4*T] ). 
Motivation for this comes, for example, from a recent claim that the fine-structure 
constant might slowly evolve with time |142j . 

5.1. The low-energy effective action 

In order to understand the cosmology of the two-brane system, we derive the low-energy 
effective action using the moduli space approximation. The moduli space approximation 
gives the low-energy limit effective action for the two-brane system, i.e. for energies 
much smaller than the brane tensions. Later in this section we will compare the moduli 
space approximation to another method. 

In the static BPS solutions described in the section 3, the brane positions can be 
chosen arbitrarily. In other words, they are moduli fields. It is expected that by putting 
some matter on the branes, these moduli fields become time- dependent, or, if the matter 
is inhomogeneously distributed, space-time dependent. Thus, the first approximation 
is to replace the brane positions with space-time dependent functions. Furthermore, in 
order to allow for the gravitational zero mode, we will replace the flat space-time metric 
r}fj, v with g^ v (x a ). We assume that the evolution of these fields is slow, which means 
that we neglect terms like (dtp) 3 when constructing the low-energy effective action. 

As already mentioned, the moduli space approximation is only a good 
approximation at energies much less than the brane tension. Thus, in the moduli space 
approximation we do not recover the quadratic term in the Friedmann equation. (See 
|143j for these issues.) We are interested in the late time effects, where the corrections 
have to be small. 

Replacing rj^ with g^(x a ) in (jlOlj) and collecting all the terms one finds from the 
five-dimensional action after an integration over y: 



with a(y) given by ()105jl . The moduli <fi and a represent the location of the two branes. 
Note that the kinetic term of the field (f) has the wrong sign. This is an artifact of the 
frame we use here. As we will see below, it is possible to go to the Einstein frame with a 
simple conformal transformation, in which the sign in front of the kinetic term is correct 
for both fields. 



[T2H].fT27j. 




(161) 



with 




(162) 
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In the following we will concentrate on the BPS system with exponential 
superpotential from section 3. Let us redefine the fields according to 

4> 2 = (l - Aka 2 <p) W , a 2 = (l - Aka 2 a) W , (163) 

with (3 = 2a ^ 1 ; and then 

= Q cosh R, a = Q sinh R. (164) 
A conformal transformation g^ u = Q 2 g fJLV leads to the Einstein frame action: 



6 -(OR) 2 



" 12a 2 (0Q) 



l + 2a 2 Q 2 

(165) 

2a 2 + l v J \ K ' 

Note that in this frame both fields have the correct sign in front of the kinetic terms. 

For a — > (i.e. the Randall-Sundrum case) the Q-field decouples. This reflects the 

fact that the bulk scalar field decouples, and the only scalar degree of freedom is the 

distance between the branes. One can read off the gravitational constant to be 

16nG = 2kK 2 5 (l + 2a 2 ). (166) 

The matter sector of the action can be found easily: if matter lives on the branes, 
it "feels" the induced metric. That is, the action has the form 

S^ = S^(^g^) and S« = £g>(* 3 , (167) 

where g^f> denotes the induced metric on each branes. In going to the Einstein frame 
one gets 

S^ = S^(^ 1 ,A 2 (Q,R)g, u ) and S® = 5^(* 2 , B 2 {Q, R) 9fll/ ), (168) 

where matter now couples explicitly to the fields via the functions A and B, which we 
will give below (neglecting derivative interactions). Note that we have assumed in ()167j) 
that the matter fields on the branes do not directly couple to the bulk scalar field. 

The theory derived with the help of the moduli space approximation has the form 
of a multiscalar-tensor theory, in which matter on both branes couples differently to 
the moduli fields. We note that methods different from the moduli-space approximation 
have been used in the literature in order to obtain the low-energy effective action 
or the resulting field equations for a two-brane system (see in particular |144j - 
|149j ). Qualitatively, the features of the resulting theories agree with the moduli-space 
approximation discussed above. We will come back to this point in section 5.4. 

In the following we will discuss observational constraints imposed on the parameters 
of the theory. 

5.2. Observational constraints 

In order to constrain the theory, it is convenient to write the moduli Lagrangian in the 
form of a non-linear sigma model with kinetic terms 

Jijdftdft , (169) 
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where i = 1, 2 labels the moduli (p 1 = Q and 2 = R. The sigma model couplings are 
here 

12a 2 1 6 

lQQ ~ TT^Q 2 "' 7M " 1+2^- (170) 

Notice the potential danger of the a — > limit, the RS model, where the coupling to Q 
becomes very small. In an ordinary Brans-Dicke theory with a single field, this would 
correspond to a vanishing Brans-Dicke parameter which is ruled out experimentally. 
Here we will see that the coupling to matter is such that this is not the case. Indeed we 
can write the action expressing the coupling to ordinary matter on our brane as 

A = ai^r 1 ' 2 ^ <t), B = a(a)r X, \4>, *), (171) 

where we have neglected the derivative interaction. 
Let us introduce the parameters 

a Q = d Q In A, a R = d R In A. (172) 

We find that (A = 4/(1 + 2a 2 )) 

A = Q-^r(coshR)^, (173) 

leading to 

a 2 A 1 Atanhi? , . 

Observations constrain the parameter 

9 = ^onoLj (175) 
to be less than 10 -5 |15Uj today. We obtain therefore a bound on 

la 2 tanh 2 R 

^ = 3TT^ + 6(l + 2a 2 )- (176) 
The bound implies that 

a < 10" 2 , R < 0.2 (177) 

today. The smallness of a indicates a strongly warped bulk geometry such as an Anti-de 
Sitter space-time. In the case a = 0, we can easily interpret the bound on R. Indeed 
in that case 

tanhi? = e - fe (*-*), (178) 

i.e. this is nothing but the exponential of the radion field measuring the distance between 
the branes. Thus we find that gravity experiments require the branes to be sufficiently 
far apart. When a ^ but small, one way of obtaining a small value of R is for the 
hidden brane to become close to the would-be singularity where a(a) = 0. 

Constraints on the parameters also arise from nucleosynthesis. Nucleosynthesis 
constrains the effective number of relativistic degrees of freedom at this epoch. In brane 
worlds the energy conservation equation implies 

pa 3 7^ const., (179) 
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resulting in a different expansion rate than that given by general relativity, giving rise 
to constraints on the parameters. One finds a < 0.1 and R < 0.4 at the time of 
nucleosynthesis. 

We would like to mention that the parameter 9 can be calculated also for matter 
on the negative tension brane. Then, following the same calculations as above, it can be 
seen that the observational constraint for 9 cannot be satisfied. Thus, if the standard 
model particles are confined on the negative tension brane, the moduli have necessarily 
to be stabilised. In the following we will assume that the standard model particles 
are confined on the positive tension brane and study the cosmological evolution of the 
moduli fields. 



5.3. Cosmological implications 

The discussion in the last subsection raises an important question: the parameter a has 
to be chosen rather small for the theory to be consistent with observations. Similarly 
the field R has to be small too. The field R is dynamical and one would like to know 
if the cosmological evolution drives the field R to small values such that it is consistent 
with the observations today. Otherwise are there natural initial conditions for the field 
Rl In the following we study the cosmological evolution of the system in order to answer 
these questions. 



5.3.1. Cosmological attractor solutions The field equations for a homogeneous and 
isotropic universe can be obtained from the action. The Friedmann equation reads 

8nG , „ T s 2a 2 - l0 1 



H A = 



(pi +P2 + V eS + W cS ) + 



3 vri ' r ~ ' c " ' ' l + 2a 2 ^ l + 2a 2 

where we have defined Q = exp 0. The field equations for R and read 



R z 



(180) 



R + 3HR 



8nG 



+ a< R(pl ~ 3 Pl) + a R'(p2 - 3p 2 ) 

1 + 2a 2 



l + 2a 2 
6 



dV eS dW eS 



OR 



dR 



+ 3/70= -8ttG 



(181) 



12a 2 



dV cS | dW cS 



+ - 3 Pi) + uf{p2 - 3p 2 ) 

The coupling parameters are given by 

2a 2 (2) 2a 2 



a 



(i) 



l + 2a 2 ' 
tanhi? 



a 



(182) 
(183) 

i + 2 « 2 ' °*' = "T+25r- (184) 

We have included matter on both branes as well as potentials V c ff and W c g. We now 
concentrate on the case where matter is only on our brane. 



.(2) 



1 + 2a 2 ' 
(tanhiTT 1 
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We will discuss the inflationary epoch in more detail below. In the radiation- 
dominated epoch the trace of the energy-momentum tensor vanishes, so that R and 
quickly become constant. The scale factor scales like a(t) oc t 1 ^ 2 . 

In the matter-dominated era the solution to these equations is given by 

/ a \ -3-2q 2 /3 / l \ 2/3-4 Q 2 /27 

Pi = Pe y—j ,a = a e ( — J (185) 



la f t \ / t \ ~ 2 / 3 

+ ln-,i? = /Jo f + i Mf) » ( 186 ) 



together with 



3 (2 e \t e J \t e 

as soon as t > t e . Note that R indeed decays. This implies that small values of R 
compatible with gravitational experiments are favoured by the cosmological evolution. 
Note, however, that the size of R in the early universe is constrained by nucleosynthesis 
as well as by the CMB anisotropies. A large discrepancy between the values of R 
during nucleosynthesis and now induces a variation of the particle masses, or equivalently 
Newton's constant, which is excluded experimentally. One can show that by putting 
matter on the negative tension brane as well, the field R evolves even faster to zero 
[126 . This behaviour is reminiscent of the attractor solution in scalar-tensor theories 

In the five-dimensional picture the fact that R is driven to small values means that 
the negative tension brane is driven towards the bulk singularity. In fact, solving the 
equations numerically for more general cases suggest that R can even be negative, which 
is, in the five-dimensional description, meaningless as the negative tension brane would 
move through the bulk singularity. Thus, in order to make any further progress, one has 
to understand the bulk singularity better|. Of course, one could simply assume that 
the negative tension brane is destroyed when it hits the singularity. A more interesting 
alternative would be if the brane is repelled instead. It was speculated that this could 
be described by some effective potential in the low-energy effective action 



5.3.2. Boundary inflation driven by an inflaton field on the brane Let us consider now 
an inflationary epoch, driven by an inflaton field embedded on the positive tension brane 
J. The action for the inflaton field x 1S > i n the Einstein frame, 

1 



S inflaton ~ / V~ 9^ x 



--A\<f>,R) (dxy-A^R)V oX r 



187) 



where we have assumed a typical chaotic inflationary potential for the inflaton field. 
Explicitly, the potential reads 

„2 



V(xA,R) = Voexp (--^—A (coshi?) 



l + 2a 



v l + 2a-y **■ (188) 

It can be shown that the field R is driven towards zero during the inflationary epoch. 

The reason is the same as in the matter era: the trace of the energy-momentum tensor 

| For a — the theory is equivalent to the Randall-Sundrum model. In this case the bulk singularity 
is shifted towards the Anti-de Sitter boundary. 
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of the inflaton field is not zero. In fact, it can be shown that R decays quickly to zero 
according to 



where c is a constant depending on the initial values of x an d and the energy scale 
Vq. In any case, R decays quickly to zero and its energy density becomes negligible. 
We assume in the following that R does not play a role in the last 60 e-folds of the 
inflationary era. We are then left with an inflationary scenario with two scalar fields. It 
is clear from eq. ()188|) that the effective energy scale 



depends on the evolution of 0. For example, for n = 2, the mass of the inflaton is not 
constant, but rather depends on the evolution of 0, since we consider the problem in the 
Einstein frame in which the Planck mass is constant, but the masses of particles vary. 
This has the effect that the background solutions are slightly modified relative to the 
case in General Relativity. 

The perturbations in this scenario were discussed in |152j for the cases of n = 2 and 
n = 4. In the case of n = 2 and a = 0.01, entropy perturbations do not play a role at 
all. The amplitude of the spectrum of entropy perturbations is around one percent that 
of the curvature perturbation. Similar conclusions hold for the case of n = 4. In figure 
1 we plot the spectral index of the curvature perturbation n R versus r T = P T /1QP R (P T 
is the tensor power spectrum and Pr is the curvature perturbation power spectrum) for 
different cases of a. As can be seen, increasing a results in an increase of the tensor 
perturbation relative to the curvature perturbation and in a smaller value of n^. The 
potential with n = 4 is already under pressure observationally, and increasing a makes 
this worse. At least for the potentials considered here, a large coupling parameter a is 
not desirable. 

5.3.3. Variation of constants In a realistic brane world model, i.e. in a model which 
is based on a fundamental theory such as superstring theory, it is expected that the 
moduli fields couple to all forms of matter embedded on the brane |153j . At the classical 
level, this automatically leads to the conclusion that gauge coupling parameters, such as 
the fine structure constant, are not fundamental parameters but depend on the vacuum 
expectation values of the moduli fields themselves. Quantum mechanically, an additional 
dependence of the couplings to the moduli fields is generated by a conformal anomaly 
when switching between the Jordan and the Einstein frame |118j . It can, however, be 
shown that one-loop corrections cancel and the variation of constants is generated only 
by directly coupling the gauge fields to the bulk scalar and having moduli-dependent 
Yukawa couplings. In the following we assume that the Yukawa couplings do not depend 
on the moduli fields. The expression for 5 a em / 'c*em is, in general, a complicated function 
of both R and <fi. However, since R rapidly decays to zero, we can neglect the variation 



R(t) oc In 



1 + exp(-4ct) 






(190) 
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Figure 1. As can be seen, increasing a results in a decrease of tir, the power law index 
of the curvature perturbation power spectrum and in an increase of tt = Pt/IGPr- 
Taken from [T55] . 



of C(em due to that of R. Then, using only the variation of 0, one can find that the 
variation of ocem as a function of z during the matter dominated epoch reads 
Sa EM 16a(3 

= -«em < ~ 2 ln(l + z). 191 

ol E m 1 + 2cr 

In this equation, (3 is a (constant) parameter, describing the coupling of photons to the 
bulk scalar field (this parameter is similar to a). 

In any case, it is a combination of a and (3 which determines the amplitude of the 
variation of a em- We can now combine two observations which constrain both a and 
(3: firstly, we can take the reported variation of cxem for granted and fit the evolution 
of cxem above to the data. Secondly, the parameter 6 defined in ([175)1 . also depends 
on (3, since the masses of baryons are directly proportional to the confinement scale 
^QCD-, which depends on the moduli fields in the Einstein frame, like all energy scales§. 
The analysis was carried out in |154j and the result is shown in figure 2 under different 
assumptions for the value of 6. 

It is interesting to note that these observations indicate that /3 ~ a, which are a 
priori independent parameters. 

§ We are considering the chiral limit here. 
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Figure 2. Constraints on the parameter a and j3, using the results from quasar 
spectra (hyperbolic curves) and different values for 9. It is interesting to note that the 
observations, taken at face value, imply that a and (3 are of the same order. Taken 
from |lfi4| . 



5.3.4- Open questions In this subsection, we have discussed some of the cosmological 
consequences of brane world moduli. It is clearly a problem for the theory that 
the negative tension brane collapses during the cosmological evolution, although this 
attractor solution makes the theory more compatible with observations. Nevertheless, 
theories of this kind suffer from the existence of singularities and it is conceivable that 
stringy effects will "smooth out" this singularity. The cosmological implications of this 
might be interesting in their own right. 

Another potentially interesting question is how Kaluza-Klein corrections affect the 
above results (see |lf)5j for an approach to include Kaluza-Klein corrections). Although 
it is unlikely that these corrections prevent the collapse of the negative tension brane, 
one can imagine that they influence the time taken for this to happen, influencing the 
dependence of R(t) on t. 



5.4- Moduli versus Projective Approaches 



So far we have considered three different approaches to the brane dynamics. The first 
one consists in solving the five dimensional equations exactly (see Section 4). This is 
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only possible in a few cases, especially when no matter is present on the branes. The 
second one is the projective approach (see Section 2 and Section 3). In that case one 
writes down the four-dimensional equations on the branes. In general one cannot find 
a closed system of equations. Indeed the projected Weyl tensor is calculable only in an 
isotropic FRW universe and reduces to the dark radiation term. In the case of a bulk 
scalar field we have seen that the second derivative of the scalar field with respect to 
the bulk coordinate is unconstrained from the point of view of the brane. Thus, the 
full dynamics of the bulk scalar field is not understood. Finally, we have presented 
the moduli space approximation in the case of a two-brane system, whereby the whole 
dynamics have been reduced to a four-dimensional description involving the low-energy 
degrees of freedom. 

In the following we will show how the projective approach and the moduli space 
approximation are linked. This is only the case in the two-brane system. In the one- 
brane system where the bulk is infinite, the calculation of the projected Weyl tensor is 
an open problem. 

We will restrict the analysis to the Randall-Sundrum case, following |145j . In that 
case the low energy degrees of freedom are the graviton and the radion. It is convenient 
to choose a gauge where the radion appears explicitly in the five dimensional metric 

ds 2 = e^dy 2 + g^(y, x^)dx^dx u (192) 

and the branes are placed at y = and y = p. In following we will consider that (f)(x^) 
is a function in four dimensions whereas g^ u (y,x^) depends on all five-dimensional 
coordinates. Decomposing the extrinsic curvature 

e-<%, = E M „ + (193) 
the junction conditions become 

^ - 7 W|o = 4 i-^K + T^v) (194) 



4 2 

on the first brane of tension Ai and carrying the energy momentum tensor T^ 1 ''^. 
Together with the Einstein equations, these junction conditions specify the dynamics of 
the system. 

As we have already mentioned several times, we will be interested in the low- 
energy limit of the brane system. This is the regime where effects relevant for the CMB 
physics take place. High energy effects would imply a modification of the early Universe 
features which might leave an imprint on the primordial fluctuations. At low energy 
on the branes, it is legitimate to use a derivative expansion as long as typical length 
scales L on the branes are much larger than the curvature of the bulk I = 1/k. Using 
dimensional analysis the Hubble rate on the brane is H = 0(L _1 ). The Friedmann 
equation when p <C A leads to p = 0(L^ 2 k^ 2 ). Using A = 0(/~ 2 kJ 2 ), this implies that 

^o((L) 2 ) (195) 




Brane World Cosmology 



41 



which is much smaller than one. Hence energy densities on the brane are much lower 
than the brane tension. Within this approximation, the bulk metric can be expanded 

g^(y,x»)=a 2 (y)\g lll ,(x' i )+g$ + ...] (196) 

where at each order the tensor gffi involves only (2n) derivatives on the brane, i.e. 
corresponds to terms of order 0((l/L) 2n ). At each order in this approximation scheme, 
the function gffl(y, x^) of y is valid throughout the bulk. Hence it encapsulates non-local 
effects due to the presence of both branes. 

In the following we will analyse the zeroth and first order terms of the derivative 
expansion. At zeroth order, we obtain the equations governing the Randall-Sundrum 
model in the absence of matter on the branes. At first order in the derivative expansion, 
the components of the extrinsic curvature tensor are given by 

^ = - \KR) + v -^{D»D v <t> - \kv 2 $) 

+ (^f - + " -fiD^) + (197) 

where covariant derivatives are taken with respect to the zeroth order metric g^(x^). 
Similarly the trace part reads 

IE iie^ v 2 p 2 ^ 
Q = ^[f + -fi D2( t> + D a ct>D a <t>) - l^D^D^] (198) 

The Einstein equations lead to the first order metric being a function of the zeroth order 
terms 

9 ® = - l ^h~ l){R ^~\ 9 ^ R) 

l 2 1 2ve^ 1 1 v 2 e 2<t} 

+ -(- - 1 - 2^ ^{D.D^ + -g^D a <j>D a <j> - y -^-(D»<t>D^ 

Zi (X LCI Zi Cl 

- l-g, u D a <pD a <P) - (1 - l)E, v (199) 

— (X 

Notice that the metric g^ and <fi are not fixed yet. The equations of motion for gravity 
and follow from the junction conditions. On the first brane, they can be cast in the 
form of Einstein equations for the induced metric g^ u 

K 2 

G» = -fT^ + E» (200) 

where is the Einstein tensor on the first brane. Of course this is nothing but 
the expected Einstein equation on the first brane. On the second brane, the induced 
metric is conformally related to the induced metric on the first one f^ v = VL 2 g^ v where 
VL = exp(— e*) and 

k 2 E^ 

Gy = -^T^ + ^ (201) 

Notice that the projected Weyl tensor appears in both Einstein equations on the branes. 
Using the explicit relation between / MJ/ and g^ one can express the Einstein tensor on 
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the second brane in terms of the Einstein tensor on the first brane and tensors depending 
on the conformal factor Q. This allows one to obtain Eg from both equations 

Eg = - ~^-{T^ + (1 - tf)T„ (2) ") - h(pPD v V - 5gD 2 ^) 

+ —^-{D^D V ^ - ^^ a fD a f)] (202) 
where ^ = 1 — Q 2 and 

wm = W^Wy (203) 

Using the tracelessness of Eg, one obtains the field equation governing the dynamics 
of ty. Moreover, now that we have determined Eg, the Einstein equations on the first 
brane lead to a closed set of equations. These equations of motion are equivalent to the 
ones obtained with the moduli approximation provided one identifies 

* ~ cosh 2 (i2) 

This shows the equivalence between the moduli space approach and the projective 
approach. The solutions of the moduli space equations of motion allow one to determine 
both Eg and the first order metric g^J. Of course, this result implies that the equations 
of motion of the moduli space approach give an accurate description of the brane 
dynamics (at lowest order in the projective approach). For other work in this direction, 
we refer to [H^-fHHj and [T56]-[T6T]. 

Having shown the relationship between the projective approach and moduli space 
approximation we now discuss cosmological perturbations and the effect the brane world 
moduli have on the CMB anisotropics. 



6. Cosmological Perturbations 

In this section we turn our attention to cosmological perturbations in brane world 
scenarios. It might be that extra dimensions have an important impact on the evolution 
(and maybe on the generation) of cosmological perturbations. Therefore, cosmological 
observations might be able to constrain extra-dimensional models, particularly with the 
advent of high precision data. 

There has been a huge number of publications related to cosmological perturbations 
in brane world scenarios, see 462"j- [T89] . In this review we just touch the topic briefly 
and describe some important results, since the details are very technical. Instead we 
refer for more details to the reviews of brane world perturbations in |19Uj and [191] . 



6.1. Cosmological Perturbations: Methods 

Generally speaking, there are two points of view one can take. The first is the point of 
view from a brane observer. In this approach one starts with the Einstein equations on 
the brane and considers perturbations of them. As we have seen, the projected Weyl 
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tensor describes the influence of the bulk. We will see below that the field equations 
on the brane do not close and hence this approach is incomplete. 

The second approach takes the point of view from the bulk. Here, one can make 
two choices of coordinates: (a) the first choice of coordinates assumes the brane to be 
at rest and straight, even in the presence of perturbations (Gaussian normal coordinate 
system). It turns out that in this coordinate system the junction conditions are easy 
to implement. However, the bulk is non-static and, more annoyingly, a coordinate 
singularity appears at finite distance from the brane |41j . One way to regulate this 
singularity is to introduce a "regulator brane", on which boundary conditions have to 
be chosen. This method was used in |192j and we will come back to this later. In 
addition to the singularity, apart from the case of a de-Sitter or Minkowski brane, the 
field equations are non-separable, (b) Because of the difficulties in the Gaussian normal 
coordinate system, one may choose coordinates in which the bulk metric take a simple 
form. Choosing a proper gauge in the bulk, one can then find the general solution to 
the bulk equations [174] . The price to be paid when choosing such a coordinate system 
is that the brane moves during the cosmological evolution and the junction conditions 
are difficult to implement. 

In two-brane systems, the situation is not easier: whatever coordinates one chooses, 
there are two junction conditions to be fulfilled. 

We discuss now the brane point of view in more detail. 

6.1.1. The brane point of view We have seen that Einstein's equations on the brane 
have the form 

G^ u = 8-kGt^ - A 4 /v + >4^tMu - E^ u (205) 

The projected Weyl tensor is tracefree: E£ = 0. Let us therefore take the point of 
view that E^ v represents a new energy-momentum tensor, representing the influence of 
the bulk gravitational field. The fluid corresponding to E^ has been named "Weyl- 
fluid". The most significant contribution from this fluid is its anisotropic stress 57Tw e yi- 
Consider the off-diagonal component of the perturbed Einstein equation. Choosing the 
longitudinal gauge, we can write the perturbed metric on the brane as 

ds 2 = -(1 - 2ij)dt 2 + (1 + 2(t))dx 2 . (206) 

Assuming that the anisotropic stress of the normal matter is zero, the off-diagonal 
component of Einstein's equation leads to 

+ tfj = -8nG57i Weyl . (207) 

Therefore, even in the absence of anisotropic stress of matter, the Weyl fluid induces 
an anisotropic stress, so that <fi + ip ^ 0. In addition, there is an entropy perturbation 
induced by the Weyl fluid. These are the two effects the bulk gravitational field has 
on cosmological perturbations. In particular, it might leave a measurable signal in the 
microwave background radiation. 
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The Sachs- Wolfe effect, relating temperature fluctuations at the last scattering 
surface to metric/matter perturbations, reads now in the brane world theory we consider 

^ = ( 6 -Pi - - STTGtfTTWeyl + ^ / da a^ Wcyl (208) 

T \ T J GR 3p cdm apl 1 J 

The first term is the usual expression for the temperature fluctuation obtained in General 
Relativity. The second term includes an entropy perturbation S*, induced by the Weyl 
fluid and is related to 5p Wey i and p = p rad + p cdm by S* = 5pwc y i/Pwe y i - <W (3(p + p)). 
The last term, however, represents an integral over the history of the anisotropic stress of 
the Weyl fluid. Because there is no equation for 7Tw e yi when considering the brane alone, 
one cannot make a prediction for the temperature anisotropy of the cosmic microwave 
background radiation. For this, the full five-dimensional problem has to be solved. 

6.2. Scalar Perturbations: CMB Anisotropics 

However, some interesting progress has been made by using low-energy approximations 
discussed in the last section. As we have seen, one considers scales much below the 
brane tension and scales much larger than the AdS curvature scale. The equation for 
the Weyl tensor reads then 

V£? F = 0. (209) 

The idea is to introduce a second brane (so-called regulator brane) in the bulk, on which 
proper boundary conditions can be chosen. The Weyl tensor is then specified by the 
interbrane distance and the energy-momentum tensor of matter on both branes. The 
evolution equation for the interbrane distance is then obtained from eq. ()209|) . 
To be more concrete, consider again the five-dimensional line element 

ds 2 = e^dy 2 + g^ v (y, x^dx^dx". (210) 

The physical distance D between the branes is given by D = exp(0)Ay, where Ay is 
the coordinate distance. Remember that we have defined ^ = 1 — exp(— 2e*). Hence, \l/ 
is completely specified by the interbrane distance. As established in the last section on 
the effective action, both the projective approach and the moduli space approximation 
lead to the conclusion that the radion couples to the trace of the energy-momentum 
tensor and therefore that the theory has much in common with a scalar-tensor theory 
(see in particular eq. (|202jl ). Thus, fluctuations of the radion satisfy a wave equation 
which is similar to that of the scalar field in a scalar-tensor theory. This then specifies 
the full evolution of the Weyl tensor and, in particular, one is able to calculate the 
effective anisotropic stress 7r,5 Weyl induced by the bulk gravitational field, so that one 
can calculate the temperature fluctuations (see eq. (|208jl ) completely. 

So far, two cases have been studied in the literature. Let us describe them briefly. 

In the first case, the second brane was regarded eulator brane in order to 

deal with boundary conditions in the bulk |192j . The interbrane distance was kept fixed, 
which can be obtained by choosing p re g.brane = —Pour.brane- Then, assuming adiabatic 
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perturbations, there is only one remaining free parameter specifying the amplitude of 
fluctuations in the dark radiation: <5Cw ey i = ^Pweyi/Prad- The results are shown in the 
figure 3. 
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Figure 3. Effects of the bulk gravitational field in brane worlds for constant radion 
and adiabatic perturbations. The quantity JCwcyi quantifies the amplitude of the 
perturbations in the Weyl fluid. C* denotes here the curvature perturbation due to 
dark radiation. Taken from |192j . 

It was found in this case that perturbations on large angular scales are affected 
most (i.e. multipole numbers smaller than 200), whereas the effect of the Weyl fluid is 
small on smaller scales. 

Let us now discuss the second case studied in the literature [193 . Here, the 
branes are free to move and there is no potential energy for the radion. The negative 
tension brane contains no matter. The background evolution was studied in the last 
section. Hence, this setup reflects a two-brane system with specific matter content 
on the two branes. For the perturbations, it was assumed that initially, fluctuations 
in the radion do not contribute to the total curvature perturbation and that they are 
subsequently sourced by fluctuations in matter. The other perturbations between matter 
and radiation are assumed to be adiabatic. The results are shown in the figure 4. 

The case of Ri n i = corresponds to the ACDM case. It can be seen from the 
normalised spectra that all scales are influenced and that the spectra for R ini > are 
damped, compared to the ACDM case. 

It is not surprising that the resulting CMB power spectra look different for the 
two cases taken from the literature. In both different behaviour of the radion 

is assumed and also the nature of fluctuations in both cases is different. Thus, the 
behaviour of the projected Weyl tensor in both cases is different. 

To summarise, although some progress has been made, it is clear that brane worlds 
will not make a definite prediction for CMB perturbations. A lot will depend on the 
dynamics of the two branes, if the radion has a potential or not, if the perturbations in 
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Figure 4. Influence of the radion on the CMB anisotropies in a two brane system 
in which the radion is allowed to move. The values in the legend denote the initial 
value of the radion in Planck units. During cosmological evolution, the radion is driven 
towards zero. The spectra are normalised to COBE. Taken from [193 . 



the radion initially play a role and whether or not there is matter on the second brane. 
However, we will be able to constrain different models for the initial conditions in the 
very early universe and, maybe, different models for inflation and the cyclic universe. 

6.3. Tensor Modes 

We have seen in section 2 that, during inflation, heavy modes decay and the zero mode is 
the only one which survives. However, are heavy modes generated during the radiation 
or matter dominated epochs? As is the case in General Relativity, the zero mode remains 
frozen outside the horizon, but new effects might appear when the massive modes re- 
enter the horizon. The only way to study this effect is to study the five-dimensional 
problem. To lowest order, one can employ a "near-brane and late-time approximation" 
and study the solution of the bulk equation near the brane [194J, 195J. In this case, 
the bulk equation is separable and can be solved explicitly. It was shown that the 
massive modes decay even outside the horizon [195 . However, nothing can be said 
about the creation of heavy modes, since this is a higher-order effect. The problem was 
attacked analytically in |196j and numerically in [197J. It was shown that massive modes 
are sourced by the initial zero mode, which re-enters the horizon and is subsequently 
damped. This effect has been studied numerically in |197j and the result is shown in 
figure 5. 

7. Brane Collisions 

The brane dynamics may lead to a collision between branes, either boundary branes 
in the Randall-Sundrum model or a floating five-brane with one of the two boundary 
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Figure 5. Damping of tensor modes in the Randall-Sundrum brane world due to 
the creation of massive modes. Here, e* = po/X at the start of integration and 7 is a 
parameter specifying the brane position. Taken from |197j . 



branes for the heterotic M-theory. This possibility has been in particular investigated 
in the context of early universe cosmology [198]- |2T3j . In all these cases, the brane 
collision is not well described at the level of a five-dimensional effective field theory. 
In particular, one may envisage that the collision may be modified by the exchange of 
open strings or open membranes when the branes get closer than a distance of the order 
of the string scale. In that case, one may have a correct description of the collision. 
Several scenarios have been considered. One appealing possibility is that the singular 
behaviour at the brane collision is resolved and a finite bounce happens, i.e. the two 
branes approach within a minimal distance and rebound before flying away in opposite 
directions. Another possibility, which has been considered in heterotic M-theory under 
the name of a "small instanton" , is that the impinging brane remains stuck to the other 
brane, implying a change in the matter content. Finally the two branes may go past 
each other, exchanging some energy in the process. The resulting configuration going 
from a big crunch to a big bang situation. This is the scenario used in the cyclic model. 

We will first discuss the cyclic model, a brane-inspired four dimensional model 
advocated to be an alternative to inflation. These models are controversial and subject 
to fierce debates. In order to illustrate the possible problems occurring during a brane 
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collision, we will consider the born again scenario. Again, we will be brief and refer to 
the reviews [214]- [2*1"""] for more details. 

The subject of brane collisions is still in its infancy, deserving further study in 
view of its cosmological relevance. In particular, the question about how perturbations 
evolve during the brane collision is an important one and has sparked a lot of activity 
[2T7]-P3K]. 

7.1. The cyclic model 

The cyclic model |2()5j is an off-shoot of the ekpyrotic model iI98] - [2~"""""| whereby the big 
bang originated from the collision of two nearly BPS branes. The primordial fluctuations 
have been argued to follow from the small ripples on the colliding branes. This scenario 
has led to heated exchanges within the cosmology community. The cyclic scenario picks 
up most of the ingredients of the ekpyrotic universe within a purely four-dimensional 
description. At its most basic level, it is described by a scalar-tensor theory involving 
one scalar degree of freedom which is inspired by the radion of the Randall-Sundrum 
model. The potential can be seen in figure 6 and can be separated into four different 
regions. There are four phases. Currently the field evolves slowly along a quintessence- 
like potential leading to the acceleration of the expansion of the universe. Later the 
field will roll down the negative part of the potential, where the universe undergoes a 
slow contracting epoch. Quantum fluctuations in the field are created during this 
phase. Then the field becomes kinetically dominated and rushes towards a singularity 
at infinity, where it bounces back to the region of acceleration. The acceleration era 
is useful in diluting the entropy created in each cycle, overcoming the problems of the 
cyclic universe envisaged by Tolman in the thirties. Now this scenario relies crucially 
on the behaviour of the universe at the singular bounce. 

The bounce occurs when the field <fr is infinite corresponding to the collision of the 
branes. The description of the physics at the bounce is an open question. In the following 
section, we will describe an attempt to understand the behaviour of the brane model 
at a singular bounce. It relies on the techniques, like the moduli space approximation, 
that we have already presented. The model was dubbed the born-again brane world 
[HE]. 

7.2. Born-again brane world 

The born-again brane world is a model of the early universe, in which two boundary 
branes collide. In doing so, the brane tensions change signs. In this model, there are no 
repeated cycles and as such it has to be distinguished from the cyclic model. 

Let us focus on the Randall-Sundrum model when the boundary branes have 
detuned tensions. We have already shown that this leads to the motion of the boundary 
branes when the tensions are increased compared to the critical value corresponding to 
the BPS static configuration. In particular this leads to the collision of the two boundary 
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Figure 6. Examples of cyclic potentials. In region (a), the potential is such that the 
field leads to the present acceleration of the expansion of the universe. In region (b), 
the universe is contracting and in region (c) the universe is driven towards the brane 
collision when the field value is — oo. Taken from |2ri4| . 



branes. More precisely the Einstein frame potential reads 

V{rj) = (Tx - 1) cosh 4 (?) + (T 2 - 1) sinh 4 (?) (211) 



,2 

where r] is a normalised field related to 

77 = -In \^-= 212 

The two branes are infinitely separated when rj — and collide when t] = oo. If T\ > 1 
and Ti + T2 < 2 the potential has a maximum. For 77 large enough, the branes are 
attracted towards each other leading to a brane collision. The branes collide when 
^ = 0. Now in the brane frame of the positive tension brane the symmetry $ — > — ^ 
exchanges positive branes and negative branes. 

In a realistic setting, our universe cannot be on the negative tension brane. Hence 
before the collision our brane was the negative tension brane which became the positive 
tension brane after the collision. Before the collision the negative tension brane frame 
can be obtained from the Einstein frame putting f^ v = g> M „/$. The field $ vanishes 
linearly at the collision leading to a divergence of the Hubble rate in the Einstein metric. 
This divergence is non-existent in the brane frame of the negative tension brane. In 
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that frame, the metric goes smoothly from the pre-collision regime (a contraction in the 
Einstein frame) to the post-collision regime (an expansion in the Einstein frame). So 
it seems that by going to the brane frame where matter couples directly to gravity, one 
can extrapolate from the pre-big bang contracting phase to the big bang epoch. 

This would justify the smooth passage used in the cyclic universe but for a hitch. 
The curvature perturbation ( possesses a logarithmic divergence at the collision. This 
divergence is present both in the Einstein and brane frames, indicating its fundamental 
nature. This indicates that the brane model is unstable at the collision. This divergence 
is at the origin of the fierce debate on cosmological perturbation in colliding brane models 
|2T7|-|235|. 

In conclusion, there is yet no known prescription to go through a cosmological 
singularity such as the one springing from the collision of two branes (for some recent 
progress, however, see |229j and p!32j. Even in the simplest setting such as the born- 
again brane world scenario, divergences show an intrinsic instability of the system. One 
may hope that by going to string theory new results on the resolution of cosmological 
singularities could be obtained. 

8. Summary 

In this review we have discussed aspects of brane world cosmology, covering both the 
simplest Randall-Sundrum model and models with a bulk scalar field, the latter being 
motivated by supersymmetry. We have seen that these models are formulated in such a 
way that their predictions can now be tested against precision data. We have discussed 
the early time cosmology when the Friedmann equation is modified and considered its 
effect on inflation, considered modifications to general relativity and how this constrains 
brane world models and discussed also the CMB predictions. More theoretical questions 
as well as speculative ideas concerning the very origin of the universe have also been 
considered. However, our review is certainly not complete. For example, we have 
been unable to cover some recent developments, such as brane cosmology in Gauss- 
Bonnet theory 236]- |239] and the cosmological consequences of induced gravity on the 
brane |240j- [2"4*H] . These topics give rise to interesting cosmological consequences; we 
refer to the original literature for more details. Similarly, several authors have studied 
brane models in dimensions higher than five (see e.g. |249j - [2*K4*j ). In this case the 
compactification is different as is the recovery of Einstein gravity. Covering this is 
beyond the scope of this review. However, we noted that in six-dimensional models 
with a bulk scalar field the problem with the self-tuning scenario is alleviated 255 . 

Brane world cosmology is in its infancy and there are many open questions. Here 
we summarise some of them. 

• In the case of the Randall-Sundrum model, the homogeneous cosmological evolution 
is well understood. However, the evolution of cosmological perturbations is still not 
fully understood. This is because the effects of the bulk gravitational field, encoded 
in the projected Weyl-tensor, on CMB physics and Large Scale Structures are 
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unknown. Whilst progress has been made in this direction, as discussed in section 
6, there is still much to do. 

• For models with bulk scalar fields, we have presented some results on the 
cosmological evolution of a homogeneous brane where we assumed that the bulk 
scalar field does not vary strongly around the brane; this needs investigating. 
Furthermore, for models with two branes, the cosmology has only been explored 
in the low energy regime, for instance with the moduli space approximation. As 
yet the cosmology of the high energy regime is an open question and can only be 
explored using techniques which go beyond the low energy effective action. 

• Both the bulk scalar field as well as the interbrane distance in two brane models 
could play an important role at least during some part of the cosmological history. 
For example one of the fields could play the role of dark energy. In that case, 
it is only natural that masses of particles, as well as other parameters, will vary. 
Preliminary work in this direction was discussed in section 5. However, more work 
needs to be done in this regard. 

• The bulk singularity seems to play an important role in a cosmological setting. 
We have seen in section 5 that the negative tension brane moves towards the 
bulk singularity. Consequently cosmology requires a detailed understanding of the 
singularity. For example, the brane might be repelled by the singularity; if this 
were the case it may have cosmological consequences. More understanding of this 
is needed. 

• Brane collisions provide an exciting development into the possible origin of the 
universe. However, much work is needed on how cosmological perturbations arise 
in this scenario and how such perturbations evolve before and after the bounce 
|217j - |2H5j . The problem is again one of the singularity. Whilst progress has been 
made there are still unanswered questions and one may well have to go beyond the 
current brane description and into the string theory to understand this singularity. 

• To date brane worlds are phenomenological models and have yet to make contact 
with an underlying theory. For example, whilst they are motivated by string theory 
the current models have not arisen as an approximation to string theory. Similarly, 
they have yet to be embedded in a realistic theory of particle physics. Much work 
is needed in this direction. 

Brane cosmology has opened up new horizons, enabling theories in extra dimensions 
to be explored and tested. Many interesting questions have arisen. With new techniques 
and high precision data cosmology will play an increasingly important role in testing 
ideas beyond the standard model of particle physics. 

Acknowledgements: We are grateful to our colleagues and friends for discussions and 
comments on brane cosmology. In particular we would like to thank Sam Webster for 
his comments on a previous version. This work was supported in part by PPARC. 



Brane World Cosmology 



52 



References 

[1] Kaluza, T. 1921 Sitzungsber.Preuss.Akad. Wiss. Berlin (Math.Phys.) Kl p 966; Klein O 1926 
Z.Phys. 37 895 

Polchinski J. 1999 String Theory, Two Volumes, Cambridge University Press 
Horava P, Witten E 1996 Nucl.Phys.B460 506; ibid Nucl.Phys.B475 94 
Witten E 1996 Nucl. Phys.B 471 135 

Lukas A, Ovrut B A, Stelle K S, Waldram D 1999 Phys.Rev.D59 086001; ibid Nucl.Phys.B 552 
246 

Akama K 1982 Lect.Notes Phys 176 267 
Rubakov V A, Shaposhnikov M E 1983 Phys.Lett.B 125 136 
Visser M 1985 Phys.Lett.B 159 22 
Squires E J 1986 Phys.Lett.B 167 286 
Gibbons G W, Wiltshire D L 1987 Nucl. Phys. B717 
Arkani-Hamed N, Dimopoulos S, Dvali G 1998 Phys.Lett.B429 263 
Antoniadis I, Arkani-Hamed N, Dimopoulos S, Dvali G 1998 Phys.Lett.B436 257 
Antoniadis I 1990 Phys.Lett.B246 377 
Randall L, Sundrum R 1999 Phys.Rev.Lett. 83 3370 
Randall L, Sundrum R 1999 Phys.Rev.Lett. 83 4690 
Goldberger W D, Wise M B 1999 Phys.Rev.Lett. 83 4922 
Binetruy P, Deffayet C, Langlois D 2000 Nucl.Phys.565 269 
Lukas A, Ovrut B A, Waldram D 1999 Phys.Rev.D 61 
Cline J M, Grojean C, Servant G 1999 Phys.Rev.Lett. 83 4245 
Csaki C, Graesser M, Kolda C, Terning J 1999 Phys.Lett.B 462 34 
Arkani-Hamed N, Dimopoulos S, Kaloper N, Sundrum R 2000 Phys. Lett.B 480 193 
Forste S, Lalak Z, Lavignac S, Nilles H-P 2000 Phys.Lett.B 481 360 
Brax Ph, van de Bruck C, Class. Quant. Grav. 20, R201 (2003) 
Lukas A, Ovrut B A, Waldram D 1998 |hep-th/9812052| 
Rubakov V A 2001 Phys.Usp.44 871 
Easson D A 2000 Int.J.Mod.Phys. A16 4803 
Wands D 2002 Class.Quant.Grav.19 3403 
Langlois D 2002 Prog. Theo. Phys. Suppl. 148 181 



Padilla A, PhD thesis 2002, University of Durham, hep-th/0210217 
Gabadadze G, |hep-ph/0308112| 
Maartens R, gr-qc/0312059 
Quevedo F 2002 Class.Quant.Grav. 19 5721 
Hoyle C D, et al 2001 Phys.Rev.Lett. 86 1418 

Flanagan E E, Tye SHE Wasserman I 2000 Phys.Rev.D 62 044039 
van de Bruck C, Dorca M, Martins C, Parry M 2000 Phys.Lett.B 495 183 
Mukohyama S 1999 Phys.Lett.B 473 241 
Ida D 2000 JHEP 0009 014 

Ichiki K, Yahiro M, Takino T, Orito M, Mathews G J 2002 Phys.Rev.D 66 043521 
Mukohyama S, Shiromizu T, Maeda K-I 2000 Phys.Rev.D 62 024028 
Kraus P 1999 JHEP 9912 011 

Binetruy P, Deffayett C, Ellwanger U, Langlois D 2000 Phys.Lett.B 477 285 
Bowcock P, Charmousis C, Gregory R 2000 Class.Quant.Grav.17 4745 



Maartens R 2001 gr-qc/0101059 



Shiromizu T, Madea K, Sasaki M 2000 Phys.Rev.D 62 024012 

Wald R 1984 General Relativity University of Chicago Press; Carroll S, Space-time and Geometry: 

Introduction to General Relativity, Addison Wesley (2003) 
Israel W 1966 Nuovo Cim B44S10 1; Erratum: ibid Nuovo Cim B48 463 



Brane World Cosmology 



53 



Maartens R 2000 Phys.Rev.D 62 084023 

Maartens R, Wands D, Bassett B A, Heard I 2000 Phys.Rev.D62 041301 
Copeland E J, Liddle A R, Lidsey J E 2001 Phys.Rev.D 64 023509 
Davis S C, Perkins W B, Davis A-C, Vernon I R 2000 Phys.Rev.D63 083518 
Langlois D, Maartens R, Sasaki M, Wands D 2001 Phys.Rev.D63 084009 
Bardeen J M, Steinhardt P J, Turner M S 1983 Phys.Rev.D28 679 
Wands D, Malik K, Lyth D H, Liddle A R (2000) Phys.Rev.D62 043527 

Liddle A R, Lyth D 2000 Cosmological Inflation and Large Scale Structure Cambridge University 
Press 

Langlois D, Maartens R, Wands D 2000 Phys.Lett.B 489 259 

Gorbunov D S, Rubakov V A, Sibiryakov S M (2001) J. High Energy Phys. JHEP 10(2001)015 
Huey G, Lidsey J E 2002 Phys.Rev.D 66 043514 
Calgagni G 2003 JCAP 03 11:009 
Calgagni G 2003, hep-ph/0312246 

Calgagni G 2004, |hep-ph/0402126| 

Seery D, Taylor A 2003, astro-ph/0309512 
Tsujikawa S, Liddle A R (2004) JCAP 0403 001 

Gordon D, Wands D, Bassett B A, Maartens R (2001) Phys.Rev.D63 023506 

Ashcroft P, van de Brack C, Davis A C 2002 Phys.Rev.D 66 121302 

Wands D, Bartolo, N, Matarrese, S, Riotto A. 2002 Phys.Rev.D 66 043520 

Emparan R, Horowitz G T, Myers R C 2000 Phys.Rev.Lett. 85 499 

Guedens R, Clancy D, Liddle A R 2002 Phys.Rev.D66 043513; 

Guedens R, Clancy D, Liddle A R 2002 Phys.Rev.D66 083509 

Carr B J, Hawking S W 1974 MNRAS 168 399 

Clancy D, Guedens R, Liddle A R 2003 Phys.Rev.D68 023507 

Majumdar A S 2003 Phys.Rev.Lett. 90 031303 

Emparan R, Garcia-Bellido J, Kaloper N 2003 JHEP 0301 079 

Majumdar A S, Mukherjee N 2004 astro-ph/0403405 



Kanti P 2004 hep-ph/0402168 



Langlois D, Sorbo L, Rodriguez-Martinez M 2002 Phys.Rev.Lett. 89 171301 

Langlois D, Sorbo L 2003 Phys.Rev.D 68 084006 

Leeper E, Maartens R, Sopuerta C F 2004 Class. Quant. Grav.21 1125 

Kogan I I, Ross G G 2000 Phys.Lett.B 485 255 

Stoica H, Tye S H H, Wasserman I 2000 Phys.Lett.B 482 205 

Davis A-C, Davis S C, Perkins W B, Vernon I R 2001 Phys.Lett.B 504, 254 

Gubser S S 2001 Phys.Rev.D63 084017 

Csaki C, Erlich J, Hollowood T J, Terning J 2001 Phys.Rev.D63 065019 
Verlinde E 2000 hep-th/0008140 
Tanaka T 2002 gr-qc/0203082 

Emparan R, Fabbri A, Kaloper N 2002 JHEP 0208 043 

Emparan R, Garcia-Bellido J, Kaloper N 2002 |hep-th/0212132| 

't Hooft G 2000 hep-th/0003004; Bousso R 2002 Rev.Mod.Phys.74, 825 

Hawking S W, Hertog T, Reall H S 2000 Phys.Rev.D62 043501 

Hawking S W, Hertog T, Reall H S 2001 Phys.Rev.D63 083504 

Nojiri S, Odinstov S D 2000 Phys.Lett.B484 119 

Anchordoqui L, Nunez C, Olsen K 2000 JHEP 0010 050 

Maeda K, Wands D 2000 Phys.Rev.D62 124009 

Mennim A, Battye R A 2001 Class. Quant. Grav. 18 2171 

Kanti P, Olive K A, Pospelov M 2000 Phys.Lett.B481 386 

Langlois D, Rodriguez-Martinez M 2001 Phys.Rev. D64 123507 

Flanagan E E, Tye S H H, Wasserman I 2001 Phys.Lett.B522 155 



Brane World Cosmology 



54 



[97] Chamblin H A, Reall H S 1999 Nucl.Phys. 562 133 

[98] Kunze K E, Vazquez-Mozo M A 2002 Phys.Rev.D65 044002 

[99] Davis S C 2002 JHEP 0203 054 

[100] Davis S C 2002 JHEP 0203 058 

[101] Kanti P, Lee S, Olive K A 2002 |hep-th/0209036| 

[102] DeWolfe O, Freedman D Z, Gubser S S, Karch A 2000 Phys.Rev.D 62 046008 

[103] Brax Ph, Davis A C 2001 Phys.Lett.B 497 289 

[104] Csaki C, Erlich J, Grojean C, Hollowood T J 2000 Nucl. Phys. B584 359 

[105] Youm D 2001 Nucl. Phys. B596 289 

[106] Brax Ph, Davis A C 2001 Phys.Lett.B 513 156 

[107] Ph. Brax, A. Falkowski and Z. Lalak, Phys. Lett. B521 (2001) 105. 

[108] Ph. Brax and N. Chatillon, JHEP 0312 (2003) 026. 

[109] Mohapatra R N, Perez-Lorenzana A, de Sousa Pires C A 2000 Phys.Rev.D62 105030 

[110] Himemoto Y, Sasaki M 2002 Phys.Rev.D 63 044015 

[111] Himemoto Y, Tanaka T, Sasaki M 2002 Phys.Rev.D 65 104020 

[112] Koyama K, Takahashi K 2003 |hep-th/0301165| 

[113] Wang B, Xue L, Zhang X, Hwang W YP 2003 |hep-th/0301072| 

[114] Brax Ph, Davis A C 2001 J. High Energy Phys. JHEP 05(2001)007 

[115] Brax Ph, van de Bruck C, Davis A C 2001 J. High Energy Phys. JHEP 10(2001)026 

[116] Chiba, T 2001 |gr-qc/01 101 18| 

[117] Uzan J-P 2003 Rev.Mod.Phys. 75, 403 

[118] Fujii Y, Maeda K-I 2003 The Scalar- Tensxor Theory of Gravitation Cambridge University Press 

[119] Kiritsis E, Kofinas G, Tetradis N, Tomaras T N, Zarikas V 2002 JHEP 0302 035 

[120] Tetradis N 2002 Phys.Lett.B569 1 

[121] Brax Ph, Langlois D, Rodriguez-Martinez M 2002 Phys.Rev.D67 104022 

[122] Garriga J, Tanaka T 2000 Phys.Rev.Lett.84 2778 

[123] Chiba T 2000 Phys.Rev.D62 021502 

[124] Brax Ph, van de Bruck C, Davis A C, Rhodes C.S. 2002 Phys.Lett.B 531 135 

[125] Brax Ph, van de Bruck C, Davis A C, Rhodes C.S. 2002 Phys.Rev.D 65 121501 

[126] Brax Ph, van de Bruck C, Davis A C, Rhodes C.S. 2003 Phys.Rev.D67 023512 

[127] Cynolter G 2002 hep-th/0209152 

[128] Csaki C, Graesser M 1, Randall L, Terning J 1999 Phys.Rev.D 62 045015 

[129] Cline J, Firouzjahi H 2001 Phys.Rev.D64 023505 

[130] Cline J, Firouzjahi H 2000 Phys.Lett.B 495, 271 

[131] Charmousis C, Gregory R, Rubakov V A 2000 Phys.Rev.D62 067505 

[132] Csaki C, Graesser M L, Kribs G D 2001 Phys.Rev.D63 065002 

[133] Gen U, Sasaki M 2001 Prog.Thcor.Phys. 105 591 

[134] Grinstein B, Noltc D R, Skiba W 2001 Phys.Rev.D63 105016 

[135] Binetruy P, Deffayet C, Langlois D 2001 Nucl.Phys.B615 219 

[136] Mazumdar A, Perez-Lorenzana A 2001 Phys.Lett.B508 340 



[137] Mazumdar A, Mohapatra R N, Perez-Lorenzana A 2003 hep-ph/0310258 
[138] Vernon I R, Davis A.-C. 2004 hep-ph/0401201 
[139] Langlois D, Sorbo L 2002 Phys.Lett.B543 155 



[140] Kofman L, Martin J, Peloso M 2004 hep-ph/0401189 

[141] Lesgourgues J, Sorbo L 2003 hep-th/0310007 

[142] Webb J K, Murphy M T, Flambaum V V, Dzuba V A, Barrow J D, Churchill C W, Prochaska 

J X, Wolfe A M 2001 Phys.Rev.Lett. 87 091301 

[143] Cline J M, Vinet J 2002 JHEP 0202, 042 

[144] Khoury J, Zhang R 2002 Phys.Rev.Lett. 89 061302 

[145] Kanno S, Soda J 2002 Phys.Rev.D66 083506 

[146] Kobayashi S, Koyama K 2002 JHEP 0212 056 



Brane World Cosmology 



55 



[147] Shiromizu T, Koyama K 2002 Phys.Rev.D67 084022 

[148] Shiromizu T, Koyama K 2003 Phys.Rev.D67 104011 

[149] Wiseman T 2002 Class.Quant.Grav. 19 3083 

[150] Damour T, Esposito-Farese G 1992 Class.Quant.Grav. 9 2093 

[151] Damour T, Nordtvedt K 1993 Phys.Rev.Lett.70 2217 

[1521 Ashcroft P R, van de Bruck C, Davis A.-C. 2003 astro-ph/0310643, to appear in Phys.Rev.D 

[153] Brax Ph, van de Bruck C, Davis A-C, Rhodes C S 2003 Astrophys. Space Sci. 283 627 

[154] Palma G A, Brax Ph, Davis A-C, van de Bruck C 2003 Phys.Rev.D68 123519 

[155] Kanno S, Soda J 2003 |hep-th/0312106| 

[156] Kanno S, Soda J 2004 Gen.Rel.Grav.36 689 

[157] Shiromizu T, Koyama K, Onda S, Torii T 2003 Phys.Rev.D 68 063506 

[158] Shiromizu T, Koyama K, Torii T 2003 Phys.Rev.D68 103513 

[159] Onda S, Shiromizu T, Koyama K, Hayakawa S 2003 |hep-th/0311262| 

[160] Yoshiguchi H, Koyama K 2004 jhep-th/0403097| 



[161] Kanno S, Soda J 2003 hep-th/0312106 

[162] Mukohymama S 2000 Phys.Rev.D 62 084015 

[163] Mukohymama S 2001 Phys.Rev.D 64 064006 

[164] Mukohymama S 2000 Class.Quant.Grav. 17 4777 

[165] Kodama H, Ishibashi A, Seto O 2000 Phys.Rev.D 62 064022 

[166] Maartens R 2000 Phys.Rev.D 62 084023 

[167] Langlois D 2000 Phys.Rev.D 62 126012 

[168] van de Bruck C, Dorca M, Brandenberger R H, Lukas A 2000 Phys.Rev.D 62 123515 

[169] Koyama K, Soda J 2000 Phys.Rev.D 62 123502 

[170] Gordon C, Maartens R 2001 Phys.Rev.D63 044022 

[171] Langlois D 2001 Phys.Rev.Lett. 86 2212 

[172] Kobayashi S, Koyama K, Soda J 2001 Phys.Lett.B501 157 

[173] Bridgman H A, Malik K A, Wands D 2002 Phys.Rev.D 65 043502 

[174] Deruelle N, Dolezel T, Katz J 2001 Phys.Rev.D 63 083513 

[175] van de Bruck C, Dorca M 2000 hep-th/0012073 



[176] Dorca M, van de Bruck C 2001 Nucl.Phys.B 605 215 

[177] Neronov A, Sachs I 2001 Phys.Lett.B 513 173 

[178] Bridgman H A, Malik K A, Wands D 2001 Phys.Rev.D 63 084012 

[179] Koyama K 2002 Phys.Rev.D66 084003 

[180] Seto O, Kodama H 2001 Phys.Rev.D 63 123506 

[181] Koyama K, Soda J 2002 Phys.Rev.D65 023514 

[182] Sago N, Himemoto Y, Sasaki M 2002 Phys.Rev.D65 024014 

[183] Leong B, Dunsby P, Challinor A, Lasenby A 2002 Phys.Rev.D65 104012 

[184] Leong B, Challinor A, Maartens R, Lasenby A 2002 Phys.Rev.D66 104010 

[185] Deffayet C 2002 Phys.Rev.D 66 103504 

[186] Giudice G F, Kolb E W, Lesgourges J, Riotto A 2002 Phys.Rev. D 66 083512 

[187] Riazuelo A, Vcrnizzi F, Steer D, Durrer R 2002 hep-th/0205220 

[188] Gumjudpai B, Maartens R, Gordon C 2003 Class.Quant.Grav. 20 3295 

[189] Ringeval C, Boe hm T, Durrer R 2003 hep-th/0307T00| 

[190] Deruelle N 2003 gr-qc/0301035 

[191] Maartens R 2004 |astro-ph7 0402485 

[192] Koyama K 2003 Phys.Rev.Lett. 91 221301 

[193] Rhodes C S, van de Bruck C, Brax Ph, Davis A.-C. 2003 Phys.Rev.D68 083511 

[194] Mennim, A 2003 PhD thesis, University of Cambridge; Mennim A, Battye R A, van de Bruck C 

2003 Astrophys. Space Sci. 283 633 

[195] Battye R A, van de Bruck C, Mennim A 2003 |hep-th /0308134 to appear in Phys.Rev.D 

[196] Easther R, Langlois D, Maartens R, Wands D 2003 JCAP 0310 014 



Brane World Cosmology 



56 



[197] Hiramatsu T, Koyama K, Taruya A 2004 Phys.Lett.B578 269 

[198] Khoury J, Ovrut B A, Steinhardt P J, Turok N 2001 Phys.Rev.D64 123522 

[199] Khoury J, Ovrut B A, Seiberg N, Steinhardt P J, Turok N 2002 Phys.Rev.D65 086007 

[200] Khoury J, Ovrut B A, Steinhardt P J, Turok N 2002 Phys.Rev.D66 046005 

[201] Enqvist K, Keski-Vakkuri E, Rasanen S 2001 Phys.Rev.D 614 388 

[202] Kallosh R, Kofman L, Linde A D 2001 Phys.Rev.D 64 123523 

[203] Neronov A 2001 JHEP 0111 007 

[204] Rasanen S 2002 Nucl.Phys. B626 183 

[205] Steinhardt P J, Turok N 2002 Science 296 1436; ibid 2002 Phys.Rev.D 65 126003 

[206] Kanno S, Sasaki M, Soda J 2003 Prog.Theor.Phys.109, 357 

[207] Ovrut B A, Pantev T, Park J 2000 JHEP 0005 045 

[208] Witten E 1996 Nucl. Phys. B460 541 

[209] Horowitz G T, Polchinski J 2002 Phys.Rev.D66 103512 

[210] Copeland E J, Gray J, Lukas A, Skinner D 2002 Phys.Rev.D66 124007 

[211] Finelli F 2003 JCAP 0310 011 

[212] Gray J, Lukas A, Probert G I 2003 |hep-th/0312111| 

[213] Antunes N, Copeland E J, Hindmarsh M, Lukas A 2003 hep-th/0310103 

[214] Rasanen S 2002 astro-ph/0208282 

[215] Khoury J 2004 |astro-ph/0401579| 

[216] Turok N, Steinhardt P J 2004 |hep-th/0403020| 

[217] Durrer R, Vernizzi F 2002 Phys.Rev.D 66 083503 

[218] Martin J, Peter P, Pinto-Neto N, Schwarz D J 2003 Phys.Rev.D67 028301 

[219] Durrer R, |hep-th/0112026| 

[220] Peter P, Pinto-Neto N 2002 Phys.Rev.D 66 063509 

[221] Brandenberger R, Finelli F 2001 JHEP0111 056 

[222] Lyth D H 2002 Phys.Lett.B524 1 

[223] Hwang J-C 2002 Phys.Rev.D 65 063514 

[224] Tsujikawa S, Brandenberger R H, Finelli F 2002 Phys.Rev.D 66 083513 

[225] Gordon C, Turok N 2003 Phys.Rev.D67 123513 

[226] Martin J, Peter P, 2003 Phys.Rev.D68 103517 

[227] Cartier C, Durrer R, Copeland E J 2003 Phys.Rev.D 67 103517 

[228] Boyle L A, Steinhardt P J, Turok N 2003 hcp-th/0307170 

[229] Tolley A J, Turok N, Steinhardt P J 2003 |hep-th/0306109| 

[230] Craps B, Ovrut B A 2004 Phys.Rev.D69 066001 

[231] Martin J, Peter P 2004 Phys.Rev.Lett.92 061301 

[232] Battefeld T J, Patil S P, Brandenberger R H 2004|hep-th/0401010| 



[233] Boyle L A, Steinhardt P J, Turok N, |hep-th /0403026 

[234] Khoury J, Steinhardt P J, Turok N 2004 Phys.Rev.Lett.92 031302 

[235] Martin J, Peter P 2004 |hep-th/0403173| 

[236] Binetruy P, Charmousis C, Davis S C, Dufaux J-F 2002 Phys.Lett.B 544 183 

[237] Germani C, Sopuerta C F 2002 Phys.Rev.Lett. 88 231101 

[238] Charmousis C, Dufaux J-F 2002 Class.Quant.Grav.19 4671 

[239] Gravanis E, Willison S 2002 |hep-th/ 0209076 

[240] Dvali G, Gabadadze G, Porrati M 2000 Phys.Lett.B485 208 

[241] Collins H, Holdom B 2000 Phys.Rcv.D62 105009 

[242] Shtanov Yu V 2000|hep^th/0005193 

[243] Deffayet C, Dvali G, Gabadadze G 2002 Phys.Rev.D65 044023 

[244] Diamandis, Georgalas, Mavromatos N, Papantonopoulos E, Pappa I 2002 Int. J. Mod. Phys. A 17 
2241 

[245] Sahni V, Shtanov Yu 2002 |astro-ph/0202346| 

[246] Maeda K I, Mizuno S, Torii T 2003 Phys.Rev.D68 024033 



Brane World Cosmology 



Multamaki, Gaztanaga E, Manera M 2003 MRNAS 344 761 
Lue A, Starkman G 2003 Phys.Rev.D67 064002 
Aghababaie Y et al 2003 JHEP 0309 037 
Schwindt, J.-M. Wetterich C 2004 Phys.Lett.B 578 409 
Graesser M L, Kile J E, Wang P 2004 hep-th/0403074 
Navarro I, Santiago J 2004 hep-th/0402204 
Lee H M, Tasinato G 2004 |hep-th/040122l| 



Battye R A, Carter B, Mennim A 2003 hep-th/0312198 



Nillcs H-P, Papazoglou A, Tasinato G 2004 Nucl.Phys.B677 405 



